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Abstract 

These proceedings summarize lectures that were delivered as part of the 2002 and 
2003 Theoretical Advanced Study Institutes in elementary particle physics (TASI) at 
the University of Colorado at Boulder. They are intended to provide a pedagogical 
introduction to cosmology aimed at advanced graduate students in particle physics and 
string theory. 
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1 Introduction 



The last decade has seen an explosive increase in both the volume and the accuracy of data 
obtained from cosmological observations. The number of techniques available to probe and 
cross-check these data has similarly proliferated in recent years. 

Theoretical cosmologists have not been slouches during this time, either. However, it is 
fair to say that we have not made comparable progress in connecting the wonderful ideas 
we have to explain the early universe to concrete fundamental physics models. One of our 
hopes in these lectures is to encourage the dialogue between cosmology, particle physics, and 
string theory that will be needed to develop such a connection. 

In this paper, we have combined material from two sets of TASI lectures (given by SMC in 
2002 and MT in 2003). We have taken the opportunity to add more detail than was originally 
presented, as well as to include some topics that were originally excluded for reasons of time. 
Our intent is to provide a concise introduction to the basics of modern cosmology as given by 
the standard "ACDM" Big-Bang model, as well as an overview of topics of current research 
interest. 

In Lecture 1 we present the fundamentals of the standard cosmology, introducing evidence 
for homogeneity and isotropy and the Priedmann-Robertson- Walker models that these make 
possible. In Lecture 2 we consider the actual state of our current universe, which leads 
naturally to a discussion of its most surprising and problematic feature: the existence of dark 
energy. In Lecture 3 we consider the implications of the cosmological solutions obtained in 
Lecture 1 for early times in the universe. In particular, we discuss thermodynamics in the 
expanding universe, finite-temperature phase transitions, and baryogenesis. Finally, Lecture 
4 contains a discussion of the problems of the standard cosmology and an introduction to 
our best-formulated approach to solving them - the inflationary universe. 

Our review is necessarily superficial, given the large number of topics relevant to modern 
cosmology. More detail can be found in several excellent textbooks [1, 2, 3, 4, 5, 6, 7]. 
Throughout the lectures we have borrowed liberally (and sometimes verbatim) from earlier 
reviews of our own [8, 9, 10, 11, 12, 13, 14, 15]. 

Our metric signature is — We use units in which h — c— 1, and define the reduced 
Planck mass by Mp = (SttG)-^/^ ~ lO^^GeV. 

2 Fundamentals of the Standard Cosmology 

2.1 Homogeneity and Isotropy: The Robertson- Walker Metric 

Cosmology as the application of general relativity (GR) to the entire universe would seem a 
hopeless endeavor were it not for a remarkable fact - the universe is spatially homogeneous 
and isotropic on the largest scales. 

"Isotropy" is the claim that the universe looks the same in all direction. Direct evidence 
comes from the smoothness of the temperature of the cosmic microwave background, as we 
will discuss later. "Homogeneity" is the claim that the universe looks the same at every 
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H X A 

Figure 2.1: Geometry of a homogeneous and isotropic space. 

point. It is harder to test directly, although some evidence comes from number counts of 
galaxies. More traditionally, we may invoke the "Copernican principle," that we do not live 
in a special place in the universe. Then it follows that, since the universe appears isotropic 
around us, it should be isotropic around every point; and a basic theorem of geometry states 
that isotropy around every point implies homogeneity. 

We may therefore approximate the universe as a spatially homogeneous and isotropic 
three-dimensional space which may expand (or, in principle, contract) as a function of time. 
The metric on such a spacetime is necessarily of the Robertson- Walker (RW) form, as we 
now demonstrate.^ 

Spatial isotropy implies spherical symmetry. Choosing a point as an origin, and using 
coordinates (r, 9, 0) around this point, the spatial line element must take the form 

da''^dr^ + f\r){de^ + sin^ed(t)'') , (1) 

where /(r) is a real function, which, if the metric is to be nonsingular at the origin, obeys 
/(r) ~ r as r — s> 0. 

Now, consider figure 2.1 in the 9 = 7i/2 plane. In this figure DH = HE = r, both DE 
and 7 arc small and HA = x. Note that the two angles labeled 7 are equal because of 
homogeneity and isotropy. Now, note that 

EF^EF'^f{2r)^^f{r)(3. (2) 

Also 

AC = 7/(r + x) = AS + = 7/(r - x) + I3f{x) . (3) 

Using (2) to eliminate /9/7, rearranging (3), dividing by 2x and taking the limit x — > 00 
yields 

= (4) 
dr 2/(r) ■ ^ ' 

^One of the authors has a sentimental attachment to the following argument, since he learned it in his 
first cosmology course [16] . 
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We must solve this subject to /(r) ~ r as r — 0. It is easy to check that if /(r) is a solution 
then f{r/a) is a solution for constant a. Also, r, sinr and sinhr are all solutions. Assuming 
analyticity and writing /(r) as a power series in r it is then easy to check that, up to scaling, 
these are the only three possible solutions. 

Therefore, the most general spacetime metric consistent with homogeneity and isotropy 

is 



ds'' = -de + a^it) [dp" + p{p) (de^ + sin^ 



where the three possibilities for /(p) are 



/(p) = {sin(p), p, sinh(p)} . 



(5) 



(6) 



This is a purely geometric fact, independent of the details of general relativity. We have used 
spherical polar coordinates (p, 6', 0), since spatial isotropy implies spherical symmetry about 
every point. The time coordinate which is the proper time as measured by a comoving 
observer (one at constant spatial coordinates) , is referred to as cosmic time, and the function 
a{t) is called the scale factor. 

There are two other useful forms for the RW metric. First, a simple change of variables 
in the radial coordinate yields 



ds^ ^ -dt'^ + a^{t) 



dr^ 



1 — kr' 



+ r'^ (dO'^ + sin^ 9d(t)'^) 



(7) 



where 

r +1 if /(p) = sin(p) 

k=\ if /(p)=p . (8) 

[ -1 if /(p) = sinh(p) 

Geometrically, k describes the curvature of the spatial sections (slices at constant cosmic 
time), k = corresponds to positively curved spatial sections (locally isometric to 3- 
spheres); A; = corresponds to local flatness, and k — —1 corresponds to negatively curved 
(locally hyperbolic) spatial sections. These are all local statements, which should be expected 
from a local theory such as GR. The global topology of the spatial sections may be that of 
the covering spaces - a 3-sphere, an infinite plane or a 3-hyperboloid - but it need not be, 
as topological identifications under freely-acting subgroups of the isometry group of each 
manifold are allowed. As a specific example, the k — Q spatial geometry could apply just as 
well to a 3-torus as to an infinite plane. 

Note that we have not chosen a normalization such that Cq = 1. We are not free to 
do this and to simultaneously normalize \k\ = 1, without including explicit factors of the 
current scale factor in the metric. In the flat case, where k = 0, we can safely choose Oq = 1. 

A second change of variables, which may be applied to either (5) or (7), is to transform 
to conformal time, r, via 
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Figure 2.2: Hubble diagrams (as replotted in [17]) showing the relationship between reces- 
sional velocities of distant galaxies and their distances. The left plot shows the original data 
of Hubble [18] (and a rather unconvincing straight-line fit through it). To reassure you, the 
right plot shows much more recent data [19], using significantly more distant galaxies (note 
difference in scale). 

Applying this to (7) yields 



where we have written a(r) = ^[^(t)] as is conventional. The conformal time does not 
measure the proper time for any particular observer, but it does simplify some calculations. 

A particularly useful quantity to define from the scale factor is the Hubble parameter 
(sometimes called the Hubble constant), given by 



The Hubble parameter relates how fast the most distant galaxies are receding from us to 
their distance from us via Hubble's law. 



This is the relationship that was discovered by Edwin Hubble, and has been verified to high 
accuracy by modern observational methods (see figure 2.2). 




(10) 



V ~ Hd. 



(12) 



7 



2.2 Dynamics: The Friedmann Equations 

As mentioned, the RW metric is a purely kinematic consequence of requiring homogeneity 
and isotropy of our spatial sections. We next turn to dynamics, in the form of differential 
equations governing the evolution of the scale factor a{t). These will come from applying 
Einstein's equation, 

R^, - ^Rg^, = 8nGT^, (13) 

to the RW metric. 

Before diving right in, it is useful to consider the types of energy-momentum tensors T^^ 
we will typically encounter in cosmology. For simplicity, and because it is consistent with 
much we have observed about the universe, it is often useful to adopt the perfect fluid form 
for the energy-momentum tensor of cosmological matter. This form is 

T^. = {P + P) Uf,U, + pfif^, , ( 14) 

where is the fluid four- velocity, p is the energy density in the rest frame of the fluid and p 
is the pressure in that same frame. The pressure is necessarily isotropic, for consistency with 
the RW metric. Similarly, fluid elements will be comoving in the cosmological rest frame, so 
that the normalized four- velocity in the coordinates of (7) will be 

(1,0,0,0). (15) 

The energy-momentum tensor thus takes the form 

/P 



T = 



P9ij 

\ / 



(16) 



where gij represents the spatial metric (including the factor of a^) . 

Armed with this simplified description for matter, we are now ready to apply Einstein's 
equation (13) to cosmology. Using (7) and (14), one obtains two equations. The first is 
known as the Friedmann equation, 

(-)=—!:«- ^. (17) 

where an overdot denotes a derivative with respect to cosmic time t and i indexes all different 
possible types of energy in the universe. This equation is a constraint equation, in the sense 
that we are not allowed to freely specify the time derivative d; it is determined in terms of 
the energy density and curvature. The second equation, which is an evolution equation, is 
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It is often useful to combine (17) and (18) to obtain the acceleration equation 



^""^EIa + SpO- (19) 



a 3 

In fact, if we know the magnitudes and evolutions of the different energy density compo- 
nents Pi, the Priedmann equation (17) is sufficient to solve for the evolution uniquely. The 
acceleration equation is conceptually useful, but rarely invoked in calculations. 

The Friedmann equation relates the rate of increase of the scale factor, as encoded by 
the Hubble parameter, to the total energy density of all matter in the universe. We may use 
the Priedmann equation to define, at any given time, a critical energy density, 

for which the spatial sections must be precisely flat {k = 0). We then define the density 
parameter 

^total = — , (21) 
Pc 

which allows us to relate the total energy density in the universe to its local geometry via 

f^totai > 1 A; = +1 

^^totai = l ^ A; = (22) 

It is often convenient to define the fractions of the critical energy density in each different 
component by 

a = - . (23) 

pc 

Energy conservation is expressed in GR by the vanishing of the covariant divergence of 
the energy-momentum tensor, 

V^T'^^ = . (24) 

Applying this to our assumptions - the RW metric (7) and perfect-fluid energy-momentum 
tensor (14) - yields a single energy-conservation equation, 

p + 3H{p + p) ^0 . (25) 

This equation is actually not independent of the Priedmann and acceleration equations, but 
is required for consistency. It implies that the expansion of the universe (as specifled by H) 
can lead to local changes in the energy density. Note that there is no notion of conservation 
of "total energy," as energy can be interchanged between matter and the spacetime geometry. 

One flnal piece of information is required before we can think about solving our cosmo- 
logical equations: how the pressure and energy density are related to each other. Within the 
fluid approximation used here, we may assume that the pressure is a single- valued function of 
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the energy density p — p{p). It is often convenient to define an equation of state parameter, 
w, by 

p — wp . (26) 

This should be thought of as the instantaneous definition of the parameter w; it need repre- 
sent the full equation of state, which would be required to calculate the behavior of fluctu- 
ations. Nevertheless, many useful cosmological matter sources do obey this relation with a 
constant value of w. For example, w = corresponds to pressureless matter, or dust - any 
collection of massive non-relativistic particles would qualify. Similarly, w = 1/3 corresponds 
to a gas of radiation, whether it be actual photons or other highly relativistic species. 

A constant w leads to a great simplification in solving our equations. In particular, 
using (25), we see that the energy density evolves with the scale factor according to 

^ ^5)^ ■ (27) 

Note that the behaviors of dust {w = 0) and radiation {w = 1/3) are consistent with what 
we would have obtained by more heuristic reasoning. Consider a fixed comoving volume of 
the universe - i.e. a volume specified by fixed values of the coordinates, from which one may 
obtain the physical volume at a given time t by multiplying by a{t)^. Given a fixed number 
of dust particles (of mass m) within this comoving volume, the energy density will then scale 
just as the physical volume, i.e. as a(t)^"^, in agreement with (27), with w = 0. 

To make a similar argument for radiation, first note that the expansion of the universe 
(the increase of a(t) with time) results in a shift to longer wavelength A, or a redshift, of 
photons propagating in this background. A photon emitted with wavelength Ag at a time tg, 
at which the scale factor is = a{te) is observed today {t — to, with scale factor oq = a{to)) 
at wavelength Ao, obeying 

^ = ^^1 + .. (28) 

Ag Og 

The redshift z is often used in place of the scale factor. Because of the redshift, the energy 
density in a fixed number of photons in a fixed comoving volume drops with the physical 
volume (as for dust) and by an extra factor of the scale factor as the expansion of the universe 
stretches the wavelengths of light. Thus, the energy density of radiation will scale as a(t)~^, 
once again in agreement with (27), with w = 1/3. 

Thus far, we have not included a cosmological constant A in the gravitational equations. 
This is because it is equivalent to treat any cosmological constant as a component of the 
energy density in the universe. In fact, adding a cosmological constant A to Einstein's 
equation is equivalent to including an energy-momentum tensor of the form 

T^ = -^9^- (29) 

This is simply a perfect fiuid with energy-momentum tensor (14) with 

A 

Pa = -pA , (30) 



10 



so that the equation-of-state parameter is 

(31) 

This imphes that the energy density is constant, 

Pa = constant . (32) 

Thus, this energy is constant throughout spacetime; we say that the cosmological constant 
is equivalent to vacuum energy. 

Similarly, it is sometimes useful to think of any nonzero spatial curvature as yet another 
component of the cosmological energy budget, obeying 

3A; 

^'^"'■^ ~ 

^'^"'^^ " 8^^' ^^^^ 

so that 

yJcuvY = -1/3 . (34) 

It is not an energy density, of course; pcurv is simply a convenient way to keep track of how 
much energy density is lacking, in comparison to a flat universe. 

2.3 Flat Universes 

It is much easier to find exact solutions to cosmological equations of motion when k = 0. 
Fortunately for us, nowadays we are able to appeal to more than mathematical simplicity to 
make this choice. Indeed, as we shall see in later lectures, modern cosmological observations, 
in particular precision measurements of the cosmic microwave background, show the universe 
today to be extremely spatially flat. 

In the case of flat spatial sections and a constant equation of state parameter w, we may 
exactly solve the Friedmann equation (27) to obtain 

ait) = ao (^- j , (35) 

where is the scale factor today, unless w = —1, in which case one obtains a{t) oc e^^. 
Applying this result to some of our favorite energy density sources yields table 1. 

Note that the matter- and radiation-dominated flat universes begin with a — 0; this is a 
singularity, known as the Big Bang. We can easily calculate the age of such a universe: 

" Jo " 3il + w)Ho ■ ^ ^ 

Unless w is close to —1, it is often useful to approximate this answer by 

to ~ H-^ . (37) 

It is for this reason that the quantity Hq^ is known as the Hubble time, and provides a useful 
estimate of the time scale for which the universe has been around. 
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Type of Energy 


p{a) 


a{t) 


Dust 




^2/3 


Radiation 




^1/2 


Cosmological Constant 


constant 





Table 1: A summary of the behaviors of the most important sources of energy density in 
cosmology. The behavior of the scale factor apphes to the case of a flat universe; the behavior 
of the energy densities is perfectly general. 

2.4 Including Curvature 

It is true that we know obscrvationally that the universe today is fiat to a high degree 
of accuracy. However, it is instructive, and useful when considering early cosmology, to 
consider how the solutions we have already identified change when curvature is included. 
Since we include this mainly for illustration we will focus on the separate cases of dust-filled 
and radiation- filled FRW models with zero cosmological constant. This calculation is an 
example of one that is made much easier by working in terms of conformal time r. 

Let us first consider models in which the energy density is dominated by matter [w = 0). 
In terms of conformal time the Einstein equations become 

k + h'^ + 2h' = , (38) 

where a prime denotes a derivative with respect to conformal time and /i(r) = a' /a. These 
equations are then easily solved for /i(r) giving 

r cot(r/2) k = 1 

h{T) = I 2/t k^O . (39) 

[ coth(T/2) A; = -1 

This then yields 

{1 — cos(t) k = 1 

rV2 k = . (40) 

cosh(r) — 1 k = —1 

One may use this to derive the connection between cosmic time and conformal time, 
which here is 

r — sin(T) k = 1 

rV6 k = . (41) 

sinh(r) — r k = —1 

Next we consider models dominated by radiation (w = 1/3). In terms of conformal time 
the Einstein equations become 

3{k + h'^) = SnGpa^ 



t{r) oc < 



k + e + 2h' = -^a^ (42) 
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Solving as we did above yields 



Mr) 



a(T) oc 



COt(T) 




1 


1/r 


k = 





coth(r) 


k = 


-1 


sin(T) 


k^ 


1 


T 


k = 





sinli(T) 


k^ 


-1 



and 



t(r) oc 



1 — cos(t) 



rV2 



k = 1 
A; = 
k^-1 



(43) 



(44) 



(45) 



cosli(r) — 1 

It is straightforward to interpret these solutions by examining the behavior of the scale 
factor a(T); the qualitative features are the same for matter- or radiation-domination. In 
both cases, the universes with positive curvature (k — +1) expand from an initial singularity 
with a = 0, and later recollapse again. The initial singularity is the Big Bang, while the final 
singularity is sometimes called the Big Crunch. The universes with zero or negative curvature 
begin at the Big Bang and expand forever. This behavior is not inevitable, however; we will 
see below how it can be altered by the presence of vacuum energy. 



2.5 Horizons 

One of the most crucial concepts to master about FRW models is the existence of horizons. 
This concept will prove useful in a variety of places in these lectures, but most importantly 
in understanding the shortcomings of what we are terming the standard cosmology. 

Suppose an emitter, e, sends a light signal to an observer, o, who is at r = 0. Setting 
6 = constant and (p = constant and working in conformal time, for such radial null rays we 
have To — T — r. In particular this means that 

To - Te = Te . (46) 

Now suppose Te is bounded below by fg; for example, fg might represent the Big Bang 
singularity. Then there exists a maximum distance to which the observer can see, known as 
the particle horizon distance, given by 

rph{ro) =To-fe . (47) 

The physical meaning of this is illustrated in figure 2.3. 

Similarly, suppose To is bounded above by fo- Then there exists a limit to spacetime 
events which can be influenced by the emitter. This limit is known as the event horizon 
distance, given by 

reh(To) = To - Te , (48) 
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Figure 2.3: Particle horizons arise when the past hght cone of an observer o terminates at a 
finite conformal time. Then there will be worldlines of other particles which do not intersect 
the past of o, meaning that they were never in causal contact. 




Figure 2.4: Event horizons arise when the future light cone of an observer o terminates at a 
finite conformal time. Then there will be worldlines of other particles which do not intersect 
the future of o, meaning that they cannot possibly infiuence each other. 
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with physical meaning illustrated in figure 2.4. 

These horizon distances may be converted to proper horizon distances at cosmic time t, 
for example 



Just as the Hubble time Hq^ provides a rough guide for the age of the universe, the Hubble 
distance cHq^ provides a rough estimate of the horizon distance in a matter- or radiation- 
dominated universe. 

2.6 Geometry, Destiny and Dark Energy 

In subsequent lectures we will use what we have learned here to extrapolate back to some of 
the earliest times in the universe. We will discuss the thermodynamics of the early universe, 
and the resulting interdependency between particle physics and cosmology. However, before 
that, we would like to explore some implications for the future of the universe. 

For a long time in cosmology, it was quite commonplace to refer to the three possible 
geometries consistent with homogeneity and isotropy as closed {k = 1), open {k = —1) and 
fiat {k = 0). There were two reasons for this. First, if one considered only the universal 
covering spaces, then a positively curved universe would be a 3-sphere, which has finite 
volume and hence is closed, while a negatively curved universe would be the hyperbolic 
3- manifold Ti.^, which has infinite volume and hence is open. 

Second, with dust and radiation as sources of energy density, universes with greater than 
the critical density would ultimately collapse, while those with less than the critical density 
would expand forever, with fiat universes lying on the border between the two. for the case 
of pure dust- filled universes this is easily seen from (40) and (44). 

As we have already mentioned, GR is a local theory, so the first of these points was never 
really valid. For example, there exist perfectly good compact hyperbolic manifolds, of finite 
volume, which are consistent with all our cosmological assumptions. However, the connection 
between geometry and destiny implied by the second point above was quite reasonable as 
long as dust and radiation were the only types of energy density relevant in the late universe. 

In recent years it has become clear that the dominant component of energy density in the 
present universe is neither dust nor radiation, but rather is dark energy. This component 
is characterized by an equation of state parameter w < —1/3. We will have a lot more to 
say about this component (including the observational evidence for it) in the next lecture, 
but for now we would just like to focus on the way in which it has completely separated our 
concepts of geometry and destiny. 

For simplicity, let's focus on what happens if the only energy density in the universe is 
a cosmological constant, with w = —1. In this case, the Friedmann equation may be solved 




(49) 
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for any value of the spatial curvature parameter k. If A > then the solutions are 



a{t) 



' cosh(y|t) 
exp Q\t) 
. sinh (^t) 



A; = +1 



A; = -1 



A; = 



(50) 



where we have encountered the A; = case earlier. It is immediately clear that, in the 
t — s> oo limit, all solutions expand exponentially, independently of the spatial curvature. In 
fact, these solutions are all exactly the same spacetime - de Sitter space - just in different 
coordinate systems. These features of de Sitter space will resurface crucially when we discuss 
inflation. However, the point here is that the universe clearly expands forever in these 
spacetimes, irrespective of the value of the spatial curvature. Note, however, that not all of 
the solutions in (50) actually cover all of de Sitter space; the /c = and k = —1 solutions 
represent coordinate patches which only cover part of the manifold. 

For completeness, let us complete the description of spaces with a cosmological constant 
by considering the case A < 0. This spacetime is called Anti-de Sitter space (AdS) and it 
should be clear from the Friedmann equation that such a spacetime can only exist in a space 
with spatial curvature A; = — 1. The corresponding solution for the scale factor is 



Once again, this solution does not cover all of AdS; for a more complete discussion, see [20]. 



In the previous lecture we set up the tools required to analyze the kinematics and dynamics 
of homogeneous and isotropic cosmologies in general relativity. In this lecture we turn to the 
actual universe in which we live, and discuss the remarkable properties cosmologists have 
discovered in the last ten years. Most remarkable among them is the fact that the universe 
is dominated by a uniformly-distributed and slowly-varying source of "dark energy," which 
may be a vacuum energy (cosmological constant), a dynamical field, or something even more 
dramatic. 

3.1 Matter: Ordinary and Dark 

In the years before we knew that dark energy was an important constituent of the universe, 
and before observations of galaxy distributions and CMB anisotropies had revolutionized the 
study of structure in the universe, observational cosmology sought to measure two numbers: 
the Hubble constant Hq and the matter density parameter ^m- Both of these quantities 
remain undeniably important, even though we have greatly broadened the scope of what we 




(51) 



3 Our Universe Today and Dark Energy 
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hope to measure. The Hubble constant is often parameterized in terms of a dimensionless 
quantity h as 

Ho = lOOh km/sec/Mpc . (52) 

After years of effort, determinations of this number seem to have zeroed in on a largely 
agreed-upon value; the Hubble Space Telescope Key Project on the extragalactic distance 
scale [21] finds 

h = 0.71 ± 0.06 , (53) 

which is consistent with other methods [22], and what we will assume henceforth. 

For years, determinations of Qu based on dynamics of galaxies and clusters have yielded 
values between approximately 0.1 and 0.4, noticeably smaller than the critical density. The 
last several years have witnessed a number of new methods being brought to bear on the 
question; here we sketch some of the most important ones. 

The traditional method to estimate the mass density of the universe is to "weigh" a cluster 
of galaxies, divide by its luminosity, and extrapolate the result to the universe as a whole. 
Although clusters are not representative samples of the universe, they are sufficiently large 
that such a procedure has a chance of working. Studies applying the virial theorem to cluster 
dynamics have typically obtained values Qm — 0.2 ±0.1 [23, 24, 25]. Although it is possible 
that the global value of M/L differs appreciably from its value in clusters, extrapolations 
from small scales do not seem to reach the critical density [26]. New techniques to weigh the 
clusters, including gravitational lensing of background galaxies [27] and temperature profiles 
of the X-ray gas [28], while not yet in perfect agreement with each other, reach essentially 
similar conclusions. 

Rather than measuring the mass relative to the luminosity density, which may be different 
inside and outside clusters, we can also measure it with respect to the baryon density [29], 
which is very likely to have the same value in clusters as elsewhere in the universe, simply 
because there is no way to segregate the baryons from the dark matter on such large scales. 
Most of the baryonic mass is in the hot intracluster gas [30] , and the fraction fg^s of total 
mass in this form can be measured either by direct observation of X-rays from the gas [31] 
or by distortions of the microwave background by scattering off hot electrons (the Sunyaev- 
Zeldovich effect) [32], typically yielding 0.1 < /gas < 0.2. Since primordial nucleosynthesis 
provides a determination of Qb ~ 0.04, these measurements imply 

^^M = ^^B//ga. = 0.3±0.1 , (54) 

consistent with the value determined from mass to light ratios. 

Another handle on the density parameter in matter comes from properties of clusters 
at high redshift. The very existence of massive clusters has been used to argue in favor of 
Qm ~ 0.2 [33], and the lack of appreciable evolution of chistcrs from high redshifts to the 
present [34, 35] provides additional evidence that < 1.0. On the other hand, a recent 
measurement of the relationship between the temperature and luminosity of X-ray clusters 
measured with the XMM-Newton sateUite [36] has been interpreted as evidence for Qm near 
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unity. This last result seems at odds with a variety of other determinations, so we should 
keep a careful watch for further developments in this kind of study. 

The story of large-scale motions is more ambiguous. The peculiar velocities of galaxies are 
sensitive to the underlying mass density, and thus to Qm, but also to the "bias" describing 
the relative amplitude of fluctuations in galaxies and mass [24, 37]. Nevertheless, recent 
advances in very large redshift surveys have led to relatively firm determinations of the mass 
density; the 2df survey, for example, finds 0.1 < flu < 0.4 [38]. 

Finally, the matter density parameter can be extracted from measurements of the power 
spectrum of density fluctuations (see for example [39]). As with the CMB, predicting the 
power spectrum requires both an assumption of the correct theory and a specification of a 
number of cosmological parameters. In simple models {e.g., with only cold dark matter and 
baryons, no massive neutrinos), the spectrum can be fit (once the amplitude is normalized) by 
a single "shape parameter", which is found to be equal to F = Quh- (For more complicated 
models see [40].) Observations then yield F ~ 0.25, or Qm ~ 0.36. For a more careful 
comparison between models and observations, see [41, 42, 43, 44]. 

Thus, we have a remarkable convergence on values for the density parameter in matter: 



As we will see below, this value is in excellent agreement with that which we would determine 
indirectly from combinations of other measurements. 

As you are undoubtedly aware, however, matter comes in different forms; the matter we 
infer from its gravitational influence need not be the same kind of ordinary matter we are 
familiar with from our experience on Earth. By "ordinary matter" we mean anything made 
from atoms and their constituents (protons, neutrons, and electrons); this would include all 
of the stars, planets, gas and dust in the universe, immediately visible or otherwise. Occa- 
sionally such matter is referred to as "baryonic matter" , where "baryons" include protons, 
neutrons, and related particles (strongly interacting particles carrying a conserved quantum 
number known as "baryon number"). Of course electrons are conceptually an important 
part of ordinary matter, but by mass they are negligible compared to protons and neutrons; 
the mass of ordinary matter comes overwhelmingly from baryons. 

Ordinary baryonic matter, it turns out, is not nearly enough to account for the observed 
matter density. Our current best estimates for the baryon density [45, 46] yield 



where these error bars are conservative by most standards. This determination comes from 
a variety of methods: direct counting of baryons (the least precise method), consistency 
with the CMB power spectrum (discussed later in this lecture), and agreement with the 
predictions of the abundances of light elements for Big-Bang nucleosynthesis (discussed in 
the next lecture). Most of the matter density must therefore be in the form of non-baryonic 
dark matter, which we will abbreviate to simply "dark matter". (Baryons can be dark, 
but it is increasingly common to reserve the terminology for the non-baryonic component.) 



0.1 <^M< 0.4 . 



(55) 



Ob = 0.04 ± 0.02 , 



(56) 
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Essentially every known particle in the Standard Model of particle physics has been ruled out 
as a candidate for this dark matter. One of the few things we know about the dark matter 
is that is must be "cold" — not only is it non-relativistic today, but it must have been that 
way for a very long time. If the dark matter were "hot" , it would have free-streamed out 
of overdense regions, suppressing the formation of galaxies. The other thing wc know about 
cold dark matter (CDM) is that it should interact very weakly with ordinary matter, so as to 
have escaped detection thus far. In the next lecture we will discuss some currently popular 
candidates for cold dark matter. 

3.2 Supernovae and the Accelerating Universe 

The great story of fin de siede cosmology was the discovery that matter does not domi- 
nate the universe; we need some form of dark energy to explain a variety of observations. 
The first direct evidence for this finding came from studies using Type la supernovae as 
"standardizable candles," which we now examine. For more detailed discussion of both the 
observational situation and the attendant theoretical problems, see [48, 49, 8, 50, 51, 15]. 

Supernovae are rare — perhaps a few per century in a Milky- Way-sized galaxy — but 
modern telescopes allow observers to probe very deeply into small regions of the sky, covering 
a very large number of galaxies in a single observing run. Supernovae are also bright, 
and Type la's in particular all seem to be of nearly uniform intrinsic luminosity (absolute 
magnitude M ~ —19.5, typically comparable to the brightness of the entire host galaxy in 
which they appear) [52]. They can therefore be detected at high redshifts {z ^ 1), allowing 
in principle a good handle on cosmological effects [53, 54]. 

The fact that all SNe la are of similar intrinsic luminosities fits well with our under- 
standing of these events as explosions which occur when a white dwarf, onto which mass is 
gradually accreting from a companion star, crosses the Chandrasekhar limit and explodes. 
(It should be noted that our understanding of supernova explosions is in a state of develop- 
ment, and theoretical models are not yet able to accurately reproduce all of the important 
features of the observed events. See [55, 56, 57] for some recent work.) The Chandrasekhar 
limit is a nearly-universal quantity, so it is not a surprise that the resulting explosions are of 
nearly-constant luminosity. However, there is still a scatter of approximately 40% in the peak 
brightness observed in nearby supernovae, which can presumably be traced to differences in 
the composition of the white dwarf atmospheres. Even if we could collect enough data that 
statistical errors could be reduced to a minimum, the existence of such an uncertainty would 
cast doubt on any attempts to study cosmology using SNe la as standard candles. 

Fortunately, the observed differences in peak luminosities of SNe la are very closely 
correlated with observed differences in the shapes of their light curves: dimmer SNe decline 
more rapidly after maximum brightness, while brighter SNe decline more slowly [58, 59, 60]. 
There is thus a one-parameter family of events, and measuring the behavior of the light curve 
along with the apparent luminosity allows us to largely correct for the intrinsic differences 
in brightness, reducing the scatter from 40% to less than 15% — sufficient precision to 
distinguish between cosmological models. (It seems likely that the single parameter can 
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Figure 3.5: Hubble diagram from the Supernova Cosmology Project, as of 2003 [70]. 

be traced to the amount of ^^Ni produced in the supernova explosion; more nickel implies 
both a higher peak luminosity and a higher temperature and thus opacity, leading to a slower 
decline. It would be an exaggeration, however, to claim that this behavior is well-understood 
theoretically.) 

Following pioneering work reported in [61], two independent groups undertook searches 
for distant supernovae in order to measure cosmological parameters: the High-Z Supernova 
Team [62, 63, 64, 65, 66], and the Supernova Cosmology Project [67, 68, 69, 70]. A plot of 
redshift vs. corrected apparent magnitude from the original SCP data is shown in Figure 3.5. 
The data are much better fit by a universe dominated by a cosmological constant than by a 
flat matter-dominated model. In fact the supernova results alone allow a substantial range 
of possible values of VLyi and f^A! however, if we think we know something about one of these 
parameters, the other will be tightly constrained. In particular, if Vtyi ~ 0.3, we obtain 

~ 0.7 . (57) 
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This corresponds to a vacuum energy density 

PA ~ 10-^ erg/cm^ ~ (10"^ eV)^ . (58) 

Thus, the supernova studies have provided direct evidence for a nonzero value for Einstein's 
cosmological constant. 

Given the significance of these results, it is natural to ask what level of confidence we 
should have in them. There are a number of potential sources of systematic error which 
have been considered by the two teams; see the original papers [63, 64, 69] for a thorough 
discussion. Most impressively, the universe implied by combining the supernova results with 
direct determinations of the matter density is spectacularly confirmed by measurements 
of the cosmic microwave background, as we discuss in the next section. Needless to say, 
however, it would be very useful to have a better understanding of both the theoretical basis 
for Type la luminosities, and experimental constraints on possible systematic errors. Future 
experiments, including a proposed satellite dedicated to supernova cosmology [71], will both 
help us improve our understanding of the physics of supernovae and allow a determination 
of the distance/redshift relation to sufficient precision to distinguish between the effects of 
a cosmological constant and those of more mundane astrophysical phenomena. 

3.3 The Cosmic Microwave Background 

Most of the radiation we observe in the universe today is in the form of an almost isotropic 
blackbody spectrum, with temperature approximately 2.7K, known as the Cosmic Microwave 
Background (CMB). The small angular fluctuations in temperature of the CMB reveal a great 
deal about the constituents of the universe, as we now discuss. 

We have mentioned several times the way in which a radiation gas evolves in and sources 
the evolution of an expanding FRW universe. It should be clear from the differing evolution 
laws for radiation and dust that as one considers earlier and earlier times in the universe, 
with smaller and smaller scale factors, the ratio of the energy density in radiation to that in 
matter grows proportionally to l/a{t). Furthermore, even particles which are now massive 
and contribute to matter used to be hotter, and at sufficiently early times were relativistic, 
and thus contributed to radiation. Therefore, the early universe was dominated by radiation. 

At early times the CMB photons were easily energetic enough to ionize hydrogen atoms 
and therefore the universe was filled with a charged plasma (and hence was opaque). This 
phase lasted until the photons redshifted enough to allow protons and electrons to combine, 
during the era of recombination. Shortly after this time, the photons decoupled from the 
now-neutral plasma and free-streamed through the universe. 

In fact, the concept of an expanding universe provides us with a clear explanation of the 
origin of the CMB. Blackbody radiation is emitted by bodies in thermal equilibrium. The 
present universe is certainly not in this state, and so without an evolving spacetime we would 
have no explanation for the origin of this radiation. However, at early times, the density and 
energy densities in the universe were high enough that matter was in approximate thermal 
equilibrium at each point in space, yielding a blackbody spectrum at early times. 
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We will have more to say about thermodynamics in the expanding universe in our next 
lecture. However, we should point out one crucial thermodynamic fact about the CMB. 
A blackbody distribution, such as that generated in the early universe, is such that at 
temperature T, the energy flux in the frequency range [u, u + du] is given by the Planck 
distribution 

P{v,T)dv = Stt/. (^)' ^,J^_^ dv , (59) 

where h is Planck's constant and k is the Boltzmann constant. Under a rescaling v — > au, 
with a^constant, the shape of the spectrum is unaltered if T — > T/a. We have already seen 
that wavelengths are stretched with the cosmic expansion, and therefore that frequencies 
will scale inversely due to the same effect. We therefore conclude that the effect of cosmic 
expansion on an initial blackbody spectrum is to retain its blackbody nature, but just at 
lower and lower temperatures, 

Tocl/a. (60) 

This is what we mean when we refer to the universe cooling as it expands. (Note that this 
strict scaling may be altered if energy is dumped into the radiation background during a 
phase transition, as we discuss in the next lecture.) 

The CMB is not a perfectly isotropic radiation bath. Deviations from isotropy at the 
level of one part in 10^ have developed over the last decade into one of our premier precision 
observational tools in cosmology. The small temperature anisotropies on the sky are usually 
analyzed by decomposing the signal into spherical harmonics via 

AT 

— = Y^aimYUeA) , (61) 

l,m 

where aim are expansion coefficients and 6 and (j) are spherical polar angles on the sky. 
Defining the power spectrum by 

Ci = {\a^'), (62) 

it is conventional to plot the quantity /(/ + l)Ci against / in a famous plot that is usually 
referred to as the CMB power spectrum. An example is shown in figure (3.6), which shows 
the measurements of the CMB anisotropy from the recent WMAP sateUite, as well as a 
theoretical model (solid line) that fits the data rather well. 

These fiuctuations in the microwave background are useful to cosmologists for many 
reasons. To understand why, we must comment briefiy on why they occur in the first place. 
Matter today in the universe is clustered into stars, galaxies, clusters and superclusters of 
galaxies. Our understanding of how large scale structure developed is that initially small 
density perturbations in our otherwise homogeneous universe grew through gravitational 
instability into the objects we observe today. Such a picture requires that from place to place 
there were small variations in the density of matter at the time that the CMB first decoupled 
from the photon-baryon plasma. Subsequent to this epoch, CMB photons propagated freely 
through the universe, nearly unaffected by anything except the cosmic expansion itself. 
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Figure 3.6: The CMB power spectrum from the WMAP satelhte [72]. The error bars on this 
plot are 1-a and the sohd hne represents the best-fit cosmological model [73] . Also shown is 
the correlation between the temperature anisotropics and the (E-mode) polarization. 
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However, at the time of their decouphng, different photons were released from regions of 
space with slightly different gravitational potentials. Since photons redshift as they climb 
out of gravitational potentials, photons from some regions redshift slightly more than those 
from other regions, giving rise to a small temperature anisotropy in the CMB observed 
today. On smaller scales, the evolution of the plasma has led to intrinsic differences in the 
temperature from point to point. In this sense the CMB carries with it a fingerprint of the 
initial conditions that ultimately gave rise to structure in the universe. 

One very important piece of data that the CMB fluctuations give us is the value of iltotai- 
Consider an overdense region of size /?, which therefore contracts under self-gravity over a 
timescale R (recall c = 1). If i? ^ -f^CMB then the region will not have had time to collapse 
over the lifetime of the universe at last scattering, li R <^ -^cmb then collapse will be well 
underway at last scattering, matter will have had time to fall into the resulting potential 
well and cause a resulting rise in temperature which, in turn, gives rise to a restoring force 
from photon pressure, which acts to damps out the inhomogeneity. 

Clearly, therefore, the maximum anisotropy will be on a scale which has had just enough 
time to collapse, but not had enough time to equilibrate - i? ~ iJj^g. This means that 
we expect to see a peak in the CMB power spectrum at an angular size corresponding to 
the horizon size at last scattering. Since we know the physical size of the horizon at last 
scattering, this provides us with a ruler on the sky. The corresponding angular scale will 
then depend on the spatial geometry of the universe. For a flat universe {k = 0, Ototai = 1) 
we expect a peak at / ~ 220 and, as can be seen in figure (3.6), this is in excellent agreement 
with observations. 

Beyond this simple heuristic description, careful analysis of all of the features of the CMB 
power spectrum (the positions and heights of each peak and trough) provide constraints on 
essentially all of the cosmological parameters. As an example we consider the results from 
WMAP [73] . For the total density of the universe they find 

0.98 < Qtotai < 1-08 (63) 

at 95% confidence - as mentioned, strong evidence for a fiat universe. Nevertheless, there 
is still some degeneracy in the parameters, and much tighter constraints on the remaining 
values can be derived by assuming either an exactly flat universe, or a reasonable value of 
the Hubble constant. When for example we assume a fiat universe, we can derive values for 
the Hubble constant, matter density (which then implies the vacuum energy density), and 
baryon density: 

h = 0.72 ±0.05 
Qm = 1 - f^A = 0.29 ± 0.07 

Qb = 0.047 ± 0.006 . 

If we instead assume that the Hubble constant is given by the value determined by the HST 
key project (53), we can derive separate tight constraints on Qm and Qa; these are shown 
graphically in Figure 3.7, along with constraints from the supernova experiments. 
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Figure 3.7: Observational constraints in the ^m-^a plane. The wide green contours represent 
constraints from supernovae, the vertical blue contours represent constraints from the 2dF 
galaxy survey, and the small orange contours represent constraints from WMAP observations 
of CMB anisotropics when a prior on the Hubble parameter is included. Courtesy of Licia 
Verde; see [74] for details. 

Taking all of the data together, we obtain a remarkably consistent picture of the current 
constituents of our universe: 





= 0.04 


DM 


= 0.26 


^^A 


= 0.7 . 



Our sense of accomplishment at having measured these numbers is substantial, although it 
is somewhat tempered by the realization that we don't understand any of them. The baryon 
density is mysterious due to the asymmetry between baryons and antibaryons; as far as dark 
matter goes, of course, we have never detected it directly and only have promising ideas as 
to what it might be. Both of these issues will be discussed in the next lecture. The biggest 
mystery is the vacuum energy; we now turn to an exploration of why it is mysterious and 
what kinds of mechanisms might be responsible for its value. 
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3.4 The Cosmological Constant Problem(s) 

In classical general relativity the cosmological constant A is a completely free parameter. It 
has dimensions of [length] (while the energy density pA has units [energy/volume]), and 
hence defines a scale, while general relativity is otherwise scale-free. Indeed, from purely 
classical considerations, we can't even say whether a specific value of A is "large" or "small" ; 
it is simply a constant of nature we should go out and determine through experiment. 

The introduction of quantum mechanics changes this story somewhat. For one thing, 
Planck's constant allows us to define the reduced Planck mass Mp ~ 10^® GeV, as well as 
the reduced Planck length 

Lp = {8nG f^ ~ 10-=^^ cm . (65) 
Hence, there is a natural expectation for the scale of the cosmological constant, namely 

^(guess) ^ ^-2 ^ ^gg^ 

or, phrased as an energy density, 

p^r'^ ~ ~ (10^^ GeV)^ ~ 10^^' erg/cm' . (67) 

We can partially justify this guess by thinking about quantum fluctuations in the vacuum. 
At all energies probed by experiment to date, the world is accurately described as a set of 
quantum fields (at higher energies it may become strings or something else). If we take 
the Fourier transform of a free quantum field, each mode of fixed wavelength behaves like a 
simple harmonic oscillator. ( "Free" means "noninteracting" ; for our purposes this is a very 
good approximation.) As we know from elementary quantum mechanics, the ground-state 
or zero-point energy of an harmonic oscillator with potential V{x) = |u;^a:^ is Eq = ^hu. 
Thus, each mode of a quantum field contributes to the vacuum energy, and the net result 
should be an integral over all of the modes. Unfortunately this integral diverges, so the 
vacuum energy appears to be infinite. However, the infinity arises from the contribution of 
modes with very small wavelengths; perhaps it was a mistake to include such modes, since 
wc don't really know what might happen at such scales. To account for our ignorance, we 
could introduce a cutoff energy, above which ignore any potential contributions, and hope 
that a more complete theory will eventually provide a physical justification for doing so. If 
this cutoff is at the Planck scale, we recover the estimate (67). 

The strategy of decomposing a free field into individual modes and assigning a zero- 
point energy to each one really only makes sense in a flat spacetime background. In curved 
spacetime we can still "renormalize" the vacuum energy, relating the classical parameter to 
the quantum value by an infinite constant. After renormalization, the vacuum energy is 
completely arbitrary, just as it was in the original classical theory. But when we use general 
relativity we are really using an effective fleld theory to describe a certain limit of quantum 
gravity. In the context of effective fleld theory, if a parameter has dimensions [mass]", we 
expect the corresponding mass parameter to be driven up to the scale at which the effective 
description breaks down. Hence, if we beheve classical general relativity up to the Planck 
scale, we would expect the vacuum energy to be given by our original guess (67). 
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However, we claim to have measured the vacuum energy (58). The observed value is 
somewhat discrepant with our theoretical estimate: 



^(obs) ^ 10-120^(guess) _ (gg) 

This is the famous 120-orders-of-magnitude discrepancy that makes the cosmological con- 
stant problem such a glaring embarrassment. Of course, it is a little unfair to emphasize the 
factor of 10^^°, which depends on the fact that energy density has units of [energy]^. We can 
express the vacuum energy in terms of a mass scale, 

Pvac = , (69) 

SO our observational result is 

Ml°^'^ - 10-3 eV . (70) 

The discrepancy is thus 

^iac'^ ~ 10-3°Af fe,"""") . (71) 

We should think of the cosmological constant problem as a discrepancy of 30 orders of 
magnitude in energy scale. 

In addition to the fact that it is very small compared to its natural value, the vacuum 
energy presents an additional puzzle: the coincidence between the observed vacuum energy 
and the current matter density. Our best-fit universe (64) features vacuum and matter 
densities of the same order of magnitude, but the ratio of these quantities changes rapidly 
as the universe expands: 

oc . (72) 



Qm Pm 

As a consequence, at early times the vacuum energy was negligible in comparison to matter 
and radiation, while at late times matter and radiation are negligible. There is only a brief 
epoch of the universe's history during which it would be possible to witness the transition 
from domination by one type of component to another. 

To date, there are not any especially promising approaches to calculating the vacuum 
energy and getting the right answer; it is nevertheless instructive to consider the example of 
supersymmetry, which relates to the cosmological constant problem in an interesting way. 
Supersymmetry posits that for each fermionic degree of freedom there is a matching bosonic 
degree of freedom, and vice- versa. By "matching" we mean, for example, that the spin-1/2 
electron must be accompanied by a spin-0 "selectron" with the same mass and charge. The 
good news is that, while bosonic fields contribute a positive vacuum energy, for fermions the 
contribution is negative. Hence, if degrees of freedom exactly match, the net vacuum energy 
sums to zero. Supersymmetry is thus an example of a theory, other than gravity, where the 
absolute zero-point of energy is a meaningful concept. (This can be traced to the fact that 
supersymmetry is a spacetime symmetry, relating particles of different spins.) 

We do not, however, live in a supersymmetric state; there is no selectron with the same 
mass and charge as an electron, or we would have noticed it long ago. If supersymmetry exists 
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in nature, it must be broken at some scale Mgusy In a theory with broken supersymmetry, 
the vacuum energy is not expected to vanish, but to be of order 

Mvac ~ M,^y , (theory) (73) 

with pvac = ^vac- What should Mgusy be? One nice feature of supersymmetry is that it helps 
us understand the hierarchy problem - why the scale of electroweak symmetry breaking is so 
much smaller than the scales of quantum gravity or grand unification. For supersymmetry 
to be relevant to the hierarchy problem, we need the supersymmetry-breaking scale to be 
just above the electroweak scale, or 

M,^y ~ 10=^ GeV . (74) 

In fact, this is very close to the experimental bound, and there is good reason to believe 
that supersymmetry will be discovered soon at Fermilab or CERN, if it is connected to 
electroweak physics. 

Unfortunately, we are left with a sizable discrepancy between theory and observation: 

Mi^'^ - . (experiment) (75) 

Compared to (71), we find that supersymmetry has, in some sense, solved the problem 
halfway (on a logarithmic scale). This is encouraging, as it at least represents a step in the 
right direction. Unfortunately, it is ultimately discouraging, since (71) was simply a guess, 
while (75) is actually a reliable result in this context; supersymmetry renders the vacuum 
energy finite and calculable, but the answer is still far away from what we need. (Subtleties in 
supergravity and string theory allow us to add a negative contribution to the vacuum energy, 
with which we could conceivably tune the answer to zero or some other small number; but 
there is no reason for this tuning to actually happen.) 

But perhaps there is something deep about supersymmetry which we don't understand, 
and our estimate Myac ~ -/^susy is simply incorrect. What if instead the correct formula were 

Mvac - (^) M.,,, ? (76) 

In other words, we are guessing that the supersymmetry-breaking scale is actually the ge- 
ometric mean of the vacuum scale and the Planck scale. Because Mp is fifteen orders of 
magnitude larger than Mgusy, and Mgusy is fifteen orders of magnitude larger than Myac, this 
guess gives us the correct answer! Unfortunately this is simply optimistic numerology; there 
is no theory that actually yields this answer (although there are speculations in this direction 
[75]). Still, the simplicity with which we can write down the formula allows us to dream that 
an improved understanding of supersymmetry might eventually yield the correct result. 

As an alternative to searching for some formula that gives the vacuum energy in terms 
of other measurable parameters, it may be that the vacuum energy is not a fundamental 
quantity, but simply our feature of our local environment. We don't turn to fundamental 
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theory for an explanation of the average temperature of the Earth's atmosphere, nor are 
we surprised that this temperature is noticeably larger than in most places in the universe; 
perhaps the cosmological constant is on the same footing. This is the idea commonly known 
as the "anthropic principle." 

To make this idea work, we need to imagine that there are many different regions of the 
universe in which the vacuum energy takes on different values; then we would expect to find 
ourselves in a region which was hospitable to our own existence. Although most humans 
don't think of the vacuum energy as playing any role in their lives, a substantially larger 
value than we presently observe would either have led to a rapid recoUapse of the universe (if 
Pvac were negative) or an inability to form galaxies (if pvac were positive). Depending on the 
distribution of possible values of p^acj one can argue that the observed value is in excellent 
agreement with what we should expect [76, 77, 78, 79, 80, 81, 82]. 

The idea of environmental selection only works under certain special circumstances, and 
we are far from understanding whether those conditions hold in our universe. In particular, 
we need to show that there can be a huge number of different domains with slightly different 
values of the vacuum energy, and that the domains can be big enough that our entire observ- 
able universe is a single domain, and that the possible variation of other physical quantities 
from domain to domain is consistent with what we observe in ours. 

Recent work in string theory has lent some support to the idea that there are a wide 
variety of possible vacuum states rather than a unique one [83, 84, 85, 86, 87, 88]. String 
theorists have been investigating novel ways to compactify extra dimensions, in which crucial 
roles are played by branes and gauge fields. By taking different combinations of extra- 
dimensional geometries, brane configurations, and gauge-field fluxes, it seems plausible that 
a wide variety of states may be constructed, with different local values of the vacuum energy 
and other physical parameters. An obstacle to understanding these purported solutions 
is the role of supersymmetry, which is an important part of string theory but needs to 
be broken to obtain a realistic universe. From the point of view of a four-dimensional 
observer, the compactiflcations that have small values of the cosmological constant would 
appear to be exactly the states alluded to earlier, where one begins with a supersymmetric 
state with a negative vacuum energy, to which supersymmetry breaking adds just the right 
amount of positive vacuum energy to give a small overall value. The necessary fine-tuning 
is accomplished simply by imagining that there are many (more than 10^°°) such states, so 
that even very unlikely things will sometimes occur. We still have a long way to go before 
we understand this possibility; in particular, it is not clear that the many states obtained 
have all the desired properties [89]. 

Even if such states are allowed, it is necessary to imagine a universe in which a large 
number of them actually exist in local regions widely separated from each other. As is 
well known, infiation works to take a small region of space and expand it to a size larger 
than the observable universe; it is not much of a stretch to imagine that a multitude of 
different domains may be separately inflated, each with different vacuum energies. Indeed, 
models of infiation generally tend to be eternal, in the sense that the universe continues to 
infiate in some regions even after infiation has ended in others [90, 91]. Thus, our observable 
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universe may be separated by inflating regions from other "universes" which have landed in 
different vacuum states; this is precisely what is needed to empower the idea of environmental 
selection. 

Nevertheless, it seems extravagant to imagine a fantastic number of separate regions of 
the universe, outside the boundary of what we can ever possibly observe, just so that wc may 
understand the value of the vacuum energy in our region. But again, this doesn't mean it 
isn't true. To decide once and for all will be extremely difficult, and will at the least require 
a much better understanding of how both string theory (or some alternative) and inflation 
operate - an understanding that we will undoubtedly require a great deal of experimental 
input to achieve. 

3.5 Dark Energy, or Worse? 

If general relativity is correct, cosmic acceleration implies there must be a dark energy density 
which diminishes relatively slowly as the universe expands. This can be seen directly from 
the Friedmann equation (17), which implies 

c? oc p + constant . (77) 

From this relation, it is clear that the only way to get acceleration (d increasing) in an 
expanding universe is if p falls off more slowly than a~^; neither matter (pm oc a~^) nor 
radiation (pr oc a~^) will do the trick. Vacuum energy is, of course, strictly constant; but 
the data are consistent with smoothly-distributed sources of dark energy that vary slowly 
with time. 

There arc good reasons to consider dynamical dark energy as an alternative to an honest 
cosmological constant. First, a dynamical energy density can be evolving slowly to zero, al- 
lowing for a solution to the cosmological constant problem which makes the ultimate vacuum 
energy vanish exactly. Second, it poses an interesting and challenging observational prob- 
lem to study the evolution of the dark energy, from which we might learn something about 
the underlying physical mechanism. Perhaps most intriguingly, allowing the dark energy to 
evolve opens the possibility of flnding a dynamical solution to the coincidence problem, if 
the dynamics are such as to trigger a recent takeover by the dark energy (independently of, 
or at least for a wide range of, the parameters in the theory). To date this hope has not 
quite been met, but dynamical mechanisms at least allow for the possibility (unlike a true 
cosmological constant). 

The simplest possibility along these lines involves the same kind of source typically in- 
voked in models of inflation in the very early universe: a scalar fleld rolling slowly in a 
potential, sometimes known as "quintessence" [92, 93, 94, 95, 96, 97]. The energy density of 
a scalar fleld is a sum of kinetic, gradient, and potential energies, 

P^ = ^<^' + ^W)' + ^(0) • (78) 
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For a homogeneous field (V0 ~ 0), the equation of motion in an expanding universe is 

+ 3i/0 + ^ = . (79) 
aq) 

If the slope of the potential V is quite flat, we will have solutions for which is nearly 
constant throughout space and only evolving very gradually with time; the energy density 
in such a configuration is 

pa constant . (80) 

Thus, a slowly-rolling scalar field is an appropriate candidate for dark energy. 

However, introducing dynamics opens up the possibility of introducing new problems, 
the form and severity of which will depend on the specific kind of model being considered. 
Most quintessence models feature scalar fields with masses of order the current Hubble 
scale, 

m<^ ~ i/o ~ 10"^^ eV . (81) 

(Fields with larger masses would typically have already rolled to the minimum of their 
potentials.) In quantum field theory, light scalar fields are unnatural; rcnormalization effects 
tend to drive scalar masses up to the scale of new physics. The well-known hierarchy problem 
of particle physics amounts to asking why the Higgs mass, thought to be of order 10^^ eV, 
should be so much smaller than the grand unification/ Planck scale, 10^^-10^^ eV. Masses of 
10~^^ eV are correspondingly harder to understand. On top of that, light scalar fields give 
rise to long-range forces and time-dependent coupling constants that should be observable 
even if couplings to ordinary matter are suppressed by the Planck scale [98, 99] ; we therefore 
need to invoke additional fine-tunings to explain why the quintessence field has not already 
been experimentally detected. 

Nevertheless, these apparent fine-tunings might be worth the price, if we were somehow 
able to explain the coincidence problem. To date, many investigations have considered scalar 
fields with potentials that asymptote gradually to zero, of the form e^^'^ or 1/0. These can 
have cosmologically interesting properties, including "tracking" behavior that makes the 
current energy density largely independent of the initial conditions [100]. They do not, 
however, provide a solution to the coincidence problem, as the era in which the scalar field 
begins to dominate is still set by finely-tuned parameters in the theory. One way to address 
the coincidence problem is to take advantage of the fact that matter/radiation equality was 
a relatively recent occurrence (at least on a logarithmic scale); if a scalar field has dynamics 
which are sensitive to the difference between matter- and radiation-dominated universes, we 
might hope that its energy density becomes constant only after matter /radiation equality. 
An approach which takes this route is fc-essence [101], which modifies the form of the kinetic 
energy for the scalar field. Instead of a conventional kinetic energy K = |(0)^, in fc-essence 
we posit a form 

K = fim^') , (82) 

where / and g are functions specified by the model. For certain choices of these functions, 
the /c-essence field naturally tracks the evolution of the total radiation energy density during 
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radiation domination, but switches to being almost constant once matter begins to dominate. 
Unfortunately, it seems necessary to choose a finely-tuned kinetic term to get the desired 
behavior [102]. 

An alternative possibility is that there is nothing special about the present era; rather, 
acceleration is just something that happens from time to time. This can be accomplished 
by oscillating dark energy [103]. In these models the potential takes the form of a decaying 
exponential (which by itself would give scaling behavior, so that the dark energy remained 
proportional to the background density) with small perturbations superimposed: 

= e-'^[l + acos(0)] . (83) 

On average, the dark energy in such a model will track that of the dominant matter/radiation 
component; however, there will be gradual oscillations from a negligible density to a dominant 
density and back, on a timescale set by the Hubble parameter, leading to occasional periods 
of acceleration. Unfortunately, in neither the A;-csscnce models nor the oscillating models do 
we have a compelling particle-physics motivation for the chosen dynamics, and in both cases 
the behavior still depends sensitively on the precise form of parameters and interactions 
chosen. Nevertheless, these theories stand as interesting attempts to address the coincidence 
problem by dynamical means. 

One of the interesting features of dynamical dark energy is that it is experimentally 
testable. In principle, different dark energy models can yield different cosmic histories, 
and, in particular, a different value for the equation of state parameter, both today and its 
redshift-dependence. Since the CMB strongly constrains the total density to be near the 
critical value, it is sensible to assume a perfectly flat universe and determine constraints on 
the matter density and dark energy equation of state; see figure (3.8) for some recent limits. 

As can be seen in (3.8), one possibility that is consistent with the data is that w < —1. 
Such a possibility violates the dominant energy condition, but possible models have been 
proposed [105]. However, such models run into serious problems when one takes them 
seriously as a particle physics theory [106, 107]. Even if one restricts one's attention to more 
conventional matter sources, making dark energy compatible with sensible particle physics 
has proven tremendously difficult. 

Given the challenge of this problem, it is worthwhile considering the possibility that 
cosmic acceleration is not due to some kind of stuff, but rather arises from new gravitational 
physics, there are a number of different approaches to this [108, 109, 110, 111, 112, 113] and 
we will not review them all here. Instead we will provide an example drawn from our own 
proposal [112]. 

As a first attempt, consider the simplest correction to the Einstein-Hilbert action. 




Here is a new parameter with units of [mass] and Cm is the Lagrangian density for matter. 

The fourth-order equations arising from this action are complicated and it is difficult to 
extract details about cosmological evolution from them. It is therefore convenient to trans- 
form from the frame used in 84, which we call the matter frame, to the Einstein frame, where 
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Figure 3.8: Constraints on the dark-energy equation-of-state parameter, as a function of 
Qu, assuming a flat universe. These hmits are derived from studies of supernovae, CMB 
anisotropies, measurements of the Hubble constant, large-scale structure, and primordial 
nucleosynthesis. From [104]. 
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the gravitational Lagrangian takes the Einstein-Hilbert form and the additional degrees of 
freedom {H and H) are represented by a fictitious scalar field (p. The details of this can 
be found in [112]. Here we just state that, performing a simultaneous redefinition of the 
time coordinate, in terms of the new metric Qf^u, our theory is that of a scalar field 
minimally coupled to Einstein gravity, and non-minimally coupled to matter, with potential 
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(85) 



Now let us first focus on vacuum cosmological solutions. The beginning of the Universe 
corresponds to i? — ^ oo and 0—^0. The initial conditions wc must specify are the initial 
values of (f) and 0', denoted as 0j and For simplicity we take 0j -C Mp. There are then 
three qualitatively distinct outcomes, depending on the value of (p'^. 

1. Eternal de Sitter. There is a critical value of (f)\ = 4>' c which just reaches the 
maximum of the potential V{4>) and comes to rest. In this case the Universe asymptotically 
evolves to a de Sitter solution. This solution requires tuning and is unstable, since any 
perturbation will induce the field to roll away from the maximum of its potential. 

2. Power-Law Acceleration. For (p\ > 4>'q, the field overshoots the maximum of V{4>) 
and the Universe evolves to late-time power-law inflation, with observational consequences 
similar to dark energy with equation-of-state parameter wde = —2/3. 

3. Future Singularity. For 0'^ < (^'q, (f) does not reach the maximum of its potential and 
rolls back down to = 0. This yields a future curvature singularity. 

In the more interesting case in which the Universe contains matter, it is possible to show 
that the three possible cosmic futures identifled in the vacuum case remain in the presence 
of matter. 

By choosing /j, ~ 10~^^ eV, the corrections to the standard cosmology only become im- 
portant at the present epoch, making this theory a candidate to explain the observed accel- 
eration of the Universe without recourse to dark energy. Since we have no particular reason 
for this choice, such a tuning appears no more attractive than the traditional choice of the 
cosmological constant. 

Clearly our choice of correction to the gravitational action can be generalized. Terms of 
the form — /i^*^"^^^/?", with n > 1, lead to similar late-time self acceleration, with behavior 
similar to a dark energy component with equation of state parameter 
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3(2n+ l)(n+ 1) 



(86) 



Clearly therefore, such modifications can easily accommodate current observational bounds [104, 
73] on the equation of state parameter —1.45 < i^de < —0.74 (95% confidence level). In the 
asymptotic regime n = 1 is ruled out at this level, while n > 2 is allowed; even n = 1 is 
permitted if we are near the top of the potential. 

Finally, any modification of the Einstein-Hilbert action must, of course, be consistent 
with the classic solar system tests of gravity theory, as well as numerous other astrophysical 
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dynamical tests. We have chosen the couphng constant n to be very small, but we have 
also introduced a new light degree of freedom. Chiba [114] has pointed out that the model 
with n = 1 is equivalent to Brans-Dicke theory with a; = in the approximation where the 
potential was neglected, and would therefore be inconsistent with experiment. It is not yet 
clear whether including the potential, or considering extensions of the original model, could 
alter this conclusion. 



4 Early Times in the Standard Cosmology 

In the first lecture we described the kinematics and dynamics of homogeneous and isotropic 
cosmologies in general relativity, while in the second we discussed the situation in our cur- 
rent universe. In this lecture we wind the clock back, using what we know of the laws of 
physics and the universe today to infer conditions in the early universe. Early times were 
characterized by very high temperatures and densities, with many particle species kept in 
(approximate) thermal equilibrium by rapid interactions. We will therefore have to move 
beyond a simple description of non-interacting "matter" and "radiation," and discuss how 
thermodynamics works in an expanding universe. 

4.1 Describing Matter 

In the first lecture we discussed how to describe matter as a perfect fiuid, described by an 
energy-momentum tensor 

T^,u ^{p + p) U^U, + pg^, , (87) 

where is the fiuid four- velocity, p is the energy density in the rest frame of the fiuid and p 
is the pressure in that same frame. The energy-momentum tensor is covariantly conserved, 

y^j^^z. = . (88) 

In a more complete description, a fiuid will be characterized by quantities in addition 
to the energy density and pressure. Many fiuids have a conserved quantity associated with 
them and so we will also introduce a number flux density N^, which is also conserved 

= . (89) 

For non-tachyonic matter N^^ is a timelike 4- vector and therefore we may decompose it as 

m = nU" . (90) 

We can also introduce an entropy flux density S^^. This quantity is not conserved, but rather 
obeys a covariant version of the second law of thermodynamics 

V^5^ > . (91) 
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Not all phenomena are successfully described in terms of such a local entropy vector (e.g., 
black holes) ; fortunately, it suffices for a wide variety of fluids relevant to cosmology. 

The conservation law for the energy-momentum tensor yields, most importantly, equa- 
tion (25), which can be thought of as the first law of thermodynamics 

dU = TdS - pdV , (92) 

with dS = 0. 

It is useful to resolve S'^ into components parallel and perpendicular to the fiuid 4- velocity 

S^' = sU^" + , (93) 

where s^f/^ = 0. The scalar s is the rest-frame entropy density which, up to an additive 
constant (that we can consistently set to zero), can be written as 

.= ^. (94) 

In addition to all these quantities, we must specify an equation of state, and we typically do 
this in such a way as to treat n and s as independent variables. 

4.2 Particles in Equilibrium 

The various particles inhabiting the early universe can be usefully characterized according to 
three criteria: in equilibrium vs. out of equilibrium (decoupled), bosonic vs. fermionic, and 
relativistic (velocities near c) vs. non-relativistic. In this section we consider species which 
are in equilibrium with the surrounding thermal bath. 

Let us begin by discussing the conditions under which a particle species will be in equilib- 
rium with the surrounding thermal plasma. A given species remains in thermal equilibrium 
as long as its interaction rate is larger than the expansion rate of the universe. Roughly 
speaking, equilibrium requires it to be possible for the products of a given reaction have 
the opportunity to recombine in the reverse reaction and if the expansion of the universe 
is rapid enough this won't happen. A particle species for which the interaction rates have 
fallen below the expansion rate of the universe is said to have frozen out or decoupled. If 
the interaction rate of some particle with the background plasma is F, it will be decoupled 
whenever 

r<^H , (95) 

where the Hubble constant H sets the cosmological timescale. 

As a good rule of thumb, the expansion rate in the early universe is "slow," and particles 
tend to be in thermal equilibrium (unless they are very weakly coupled). This can be seen 
from the Friedmann equation when the energy density is dominated by a plasma with p ~ T^; 
we then have 

H-(^)t. (96) 
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Thus, the Hubble parameter is suppressed with respect to the temperature by a factor of 
T/Mp. At extremely early times (near the Planck era, for example), the universe may 
be expanding so quickly that no species are in equilibrium; as the expansion rate slows, 
equilibrium becomes possible. However, the interaction rate F for a particle with cross- 
section cr is typically of the form 

r = n{av) , (97) 

where n is the number density and v a typical particle velocity. Since n oc a~^, the density 
of particles will eventually dip so low that equilibrium can once again no longer be main- 
tained. In our current universe, no species are in equilibrium with the background plasma 
(represented by the CMB photons). 

Now let us focus on particles in equilibrium. For a gas of wcakly-intcracting particles, we 
can describe the state in terms of a distribution function /(p), where the three-momentum 
p satisfies 

E^{p) = + \pf . (98) 

The distribution function characterizes the density of particles in a given momentum bin. 
(In general it will also be a function of the spatial position x, but we suppress that here.) 
The number density, energy density, and pressure of some species labeled i are given by 

where gi is the number of spin states of the particles. For massless photons we have = 2, 
while for a massive vector boson such as the Z we have gz — ^- In the usual accounting, 
particles and antiparticles are treated as separate species; thus, for spin-1/2 electrons and 
positrons we have g,.- = g^+ =2. In thermal equilibrium at a temperature T the particles 
will be in either Fermi-Dirac or Bose-Einstein distributions, 

/(P) = ^555^ • 

where the plus sign is for fcrmions and the minus sign for bosons. 

We can do the integrals over the distribution functions in two opposite limits: particles 
which are highly relativistic (T 3> m) or highly non-relativistic (T <^ m). The results are 
shown in table 2, in which C, is the Riemann zeta function, and C(3) ~ 1.202. 

From this table we can extract several pieces of relevant information. Relativistic parti- 
cles, whether bosons or fermions, remain in approximately equal abundances in equilibrium. 
Once they become non-relativistic, however, their abundance plummets, and becomes expo- 
nentially suppressed with respect to the relativistic species. This is simply because it becomes 
progressively harder for massive particle- antiparticle pairs to be produced in a plasma with 

r< m. 
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Table 2: Number density, energy density, and pressure, for species in thermal equilibrium. 

It is interesting to note that, although matter is much more dominant than radiation 
in the universe today, since their energy densities scale differently the early universe was 
radiation-dominated. We can write the ratio of the density parameters in matter and radi- 
ation as 

,1 + ^) ■ (101) 



^Ro Vao/ ^RO 
The redshift of matter-radiation equality is thus 

1 + ^eq = ^ 3 X 10^ . (102) 

This expression assumes that the particles that are non-relativistic today were also non- 
relativistic at Zeq, this should be a safe assumption, with the possible exception of massive 
neutrinos, which make a minority contribution to the total density. 

As we mentioned in our discussion of the CMB in the previous lecture, even decoupled 
photons maintain a thermal distribution; this is not because they are in equilibrium, but 
simply because the distribution function redshifts into a similar distribution with a lower 
temperature proportional to 1/a. We can therefore speak of the "effective temperature" of 
a relativistic species that freezes out at a temperature Tf and scale factor a/: 

Tr\a) = Tf (^) . (103) 

For example, neutrinos decouple at a temperature around 1 McV; shortly thereafter, electrons 
and positrons annihilate into photons, dumping energy (and entropy) into the plasma but 
leaving the neutrinos unaffected. Consequently, we expect a neutrino background in the 
current universe with a temperature of approximately 2K, while the photon temperature is 
3K. 

A similar effect occurs for particles which are non-relativistic at decoupling, with one 
important difference. For non-relativistic particles the temperature is proportional to the 
kinetic energy |mv^, which redshifts as We therefore have 
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;a) = r,(^)'. (104) 
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In either case we are imagining that the species freezes out while relativistic/non-relativistic 
and stays that way afterward; if it freezes out while relativistic and subsequently becomes 
non-relativistic, the distribution function will be distorted away from a thermal spectrum. 

The notion of an effective temperature allows us to define a corresponding notion of 
an effective number of relativistic degrees of freedom, which in turn permits a compact 
expression for the total relativistic energy density. The effective number of relativistic degrees 
of freedom (as far as energy is concerned) can be defined as 

s.= i:«.(f)^5^E «(f)'. (105) 

bosons lermions 

(The temperature T is the actual temperature of the background plasma, assumed to be in 
equilibrium.) Then the total energy density in all relativistic species comes from adding the 
contributions of each species, to obtain the simple formula 

P^'^g^T'. (106) 

We can do the same thing for the entropy density. From (94), the entropy density in rela- 
tivistic particles goes as rather than T^, so we define the effective number of relativistic 
degrees of freedom for entropy as 

bosons fermions 

The entropy density in relativistic species is then 

s = '%9*sT^ ■ (108) 

Numerically, g* and g^S will typically be very close to each other. In the Standard Model, 
we have 

r 100 T > 300 MeV 
g* ~ f*s < 10 300 MeV >T>1 MeV (109) 
I 3 T <l MeV . 

The events that change the effective number of relativistic degrees of freedom are the QCD 
phase transition at 300 MeV, and the annihilation of electron/positron pairs at 1 MeV. 

Because of the release of energy into the background plasma when species annihilate, it is 
only an approximation to say that the temperature goes as T oc 1/a. A better approximation 
is to say that the comoving entropy density is conserved, 

s oc . (110) 

This will hold under all forms of adiabatic evolution; entropy will only be produced at a 
process like a first-order phase transition or an out-of-equilibrium decay. (In fact, we expect 
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that the entropy production from such processes is very small compared to the total entropy, 
and adiabatic evolution is an excellent approximation for almost the entire early universe. 
One exception is inflation, discussed in the next lecture.) Combining entropy conservation 
with the expression (108) for the entropy density in relativistic species, we obtain a better 
expression for the evolution of the temperature, 

T^g:^'\-^. (Ill) 

The temperature will consistently decrease under adiabatic evolution in an expanding uni- 
verse, but it decreases more slowly when the effective number of relativistic degrees of freedom 
is diminished. 

4.3 Thermal Relics 

As we have mentioned, particles typically do not stay in equilibrium forever; eventually 
the density becomes so low that interactions become infrequent, and the particles freeze 
out. Since essentially all of the particles in our current universe fall into this category, it is 
important to study the relic abundance of decoupled species. (Of course it is also possible 
to obtain a significant relic abundance for particles which were never in thermal equilibrium; 
examples might include baryons produced by GUT baryogenesis, or axions produced by 
vacuum misalignment.) In this section we will typically neglect factors of order unity. 

We have seen that relativistic, or hot^ particles have a number density that is proportional 
to in equilibrium. Thus, a species X that freezes out while still relativistic will have a 
number density at freeze-out T/ given by 

nx{Tf) ^ T| . (112) 

Since this is comparable to the number density of photons at that time, and after freeze-out 
both photons and our species X just have their number densities dilute by a factor a{t)~^ 
as the universe expands, it is simple to see that the abundance of X particles today should 
be comparable to the abundance of CMB photons, 

nxQ ~ n^Q ~ 10^ cm~^ . (113) 

We express this number as 10^ rather than 411 since the roughness of our estimate does not 
warrant such misleading precision. The leading correction to this value is typically due to 
the production of additional photons subsequent to the decoupling of X; in the Standard 
Model, the number density of photons increases by a factor of approximately 100 between 
the electroweak phase transition and today, and a species which decouples during this period 
will be diluted by a factor of between 1 and 100 depending on precisely when it freezes out. 
So, for example, neutrinos which are light {nii, < MeV) have a number density today of 
Hi, = 115 cm~^ per species, and a corresponding contribution to the density parameter (if 
they are nevertheless heavy enough to be nonrelativistic today) of 
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(In this final expression we have secretly taken account of the missing numerical factors, 
so this is a reliable answer.) Thus, a neutrino with rrii, ~ 10~^ eV would contribute fl^, ~ 
2 X 10"^. This is large enough to be interesting without being large enough to make neutrinos 
be the dark matter. That's good news, since the large velocities of neutrinos make them 
free-stream out of overdense regions, diminishing primordial perturbations and leaving us 
with a universe which has much less structure on small scales than we actually observe. 

Now consider instead a species X which is nonrelativistic or cold at the time of decoupling. 
It is much harder to accurately calculate the relic abundance of a cold relic than a hot one, 
simply because the equilibrium abundance of a nonrelativistic species is changing rapidly with 
respect to the background plasma, and we have to be quite precise following the freeze-out 
process to obtain a reliable answer. The accurate calculation typically involves numerical 
integration of the Boltzmann equation for a network of interacting particle species; here, 
we cut to the chase and simply provide a reasonable approximate expression. If (Tq is the 
annihilation cross-section of the species X at a temperature T — mx, the final number 
density in terms of the photon density works out to be 

nx{T < TA ^n-n^ . (115) 

Since the particles are nonrelativistic when they decouple, they will certainly be nonrela- 
tivistic today, and their energy density is 

Px = mxnx . (116) 

We can plug in numbers for the Hubble parameter and photon density to obtain the density 
parameter, 

^x-^- . (117) 
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Numerically, when h = c = 1 we have 1 GeV~ 2 x 10 cm, so the photon density today is 
~ 100 cm-3 ~ 10-=^^ GeV-=^. The Hubble constant is Hq ~ 10"^^ Q^y, and the Planck 
mass is Mp ~ 10^^ GeV, so we obtain 

It is interesting to note that this final expression is independent of the mass nix of our relic, 
and only depends on the annihilation cross-section; that's because more massive particles will 
have a lower relic abundance. Of course, this depends on how we choose to characterize our 
theory; we may use variables in which ctq is a function of mx, in which case it is reasonable 
to say that the density parameter does depend on the mass. 

The designation cold may ring a bell with many of you, for you will have heard it used in a 
cosmological context applied to cold dark matter (CDM). Let us see briefly why this is. One 
candidate for CDM is a Weakly Interacting Massive Particle (WIMP). The annihilation cross- 
section of these particles, since they are weakly interacting, should be uo ~ a^^Gp, where aw 
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is the weak coupling constant and Gp is the the Fermi constant. Using Gp ~ (300GeV) ^ 
and aw ~ 10~^, we get 

(To ~ alrGp ~ 10"^ GeV"^ . (119) 
Thus, the density parameter in such particles would be 

~ 1 ■ (120) 

In other words, a stable particle with a weak interaction cross section naturally produces a 
relic density of order the critical density today, and so provides a perfect candidate for cold 
dark matter. A paradigmatic example is provided by the lightest supersymmetric partner 
(LSP), if it is stable and supersymmetry is broken at the weak scale. Such a possibility is 
of great interest to both particle physicists and cosmologists, since it may be possible to 
produce and detect such particles in colliders and to directly detect a WIMP background 
in cryogenic detectors in underground laboratories; this will be a major experimental effort 
over the next few years [13]. 



4.4 Vacuum displacement 

Another important possibility is the existence of relics which were never in thermal equi- 
librium. An example of these will be discussed later in this lecture: the production of 
topological defects at phase transitions. Let's discuss another kind of non-thermal relic, 
which derives from what we might call "vacuum displacement." Consider the action for a 
real scalar field in curved spacetime: 



-^g^'^d,<f>d,<f>-V{<f>) 



(121) 



If we assume that is spatially homogeneous {di(j) — 0), its equation of motion in the 
Robertson- Walker metric (5) will be 

(f> + 3H<j) + V'{(l)) ^0 , (122) 

where an overdot indicates a partial derivative with respect to time, and a prime indicates 
a derivative with respect to (j). For a free massive scalar field, V{(p) = |m^0^, and (122) 
describes a harmonic oscillator with a time-dependent damping term. For H > rrirj, the field 
will be overdamped, and stay essentially constant at whatever point in the potential it finds 
itself. So let us imagine that at some time in the very early universe (when H was large) we 
had such an overdamped homogeneous scalar field, stuck at a value (p = (p^,; the total energy 
density in the field is simply the potential energy \m\(j)1. The Hubble parameter H will 



decrease to approximately when the temperature reaches T* = ^m^jyMp^ after which the 
field will be able to evolve and will begin to oscillate in its potential. The vacuum energy 
is converted to a combination of vacuum and kinetic energy which will redshift like matter. 
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as p0 oc a ^; in a particle interpretation, the field is a Bose condensate of zero- momentum 
particles. We will therefore have 



<^)-l<<^l[^) ' (123) 
which leads to a density parameter today 

A classic example of a non-thermal relic produced by vacuum displacement is the QCD 
axion, which has a typical primordial value (0) ~ /pq and a mass ~ ^qcd//pQ' where 
/pQ is the Peccei-Quinn symmetry- breaking scale and Aqcd ~ 0.3 GeV is the QCD scale [1]. 
In this case, plugging in numbers reveals 

The Peccci-Quinn scale is essentially a free parameter from a theoretical point of view, 
but experiments and astrophysical constraints have ruled out most values except for a small 
window around /pq ~ 10^^ GeV. The axion therefore remains a viable dark matter candidate 
[115, 116]. Note that, even though dark matter axions are very light (Aqcd//pq 10~^ eV), 
they are extremely non-relativistic, which can be traced to the non-thermal nature of their 
production process. (Another important way to produce axions is through the decay of axion 
cosmic strings [1, 117].) 



4.5 Primordial Nucleosynthesis 

Given our time constraints, even some of the more important concepts in cosmology cannot 
be dealt with in significant detail. We have chosen just to give a cursory treatment to 
primordial nucleosynthesis, although its importance as a crucial piece of evidence in favor of 
the big bang model and its usefulness in bounding any new physics of cosmological relevance 
cannot be overstated. 

Observations of primordial nebulae reveal abundances of the light elements unexplained 
by stellar nucleosynthesis. Although it does a great disservice to the analytic and numerical 
work required, not to mention the difficulties of measuring the abundances, we will just state 
that the study of nuclear processes in the background of an expanding cooling universe yields 
a remarkable concordance between theory and experiment. 

At temperatures below 1 MeV, the weak interactions are frozen out and neutrons and 
protons cease to intcrconvcrt. The equilibrium abundance of neutrons at this temperature 
is about 1/6 the abundance of protons (due to the slightly larger neutron mass). The 
neutrons have a finite lifetime (t„ = 890 sec) that is somewhat larger than the age of the 
universe at this epoch, t{l MeV) 1 sec, but they begin to gradually decay into protons and 
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leptons. Soon thereafter, however, we reach a temperature somewhat below 100 keV, and 
Big-Bang nucleosynthesis (BBN) begins. (The nuclear binding energy per nucleon is typically 
of order 1 MeV, so you might expect that nucleosynthesis would occur earlier; however, the 
large number of photons per nucleon prevents nucleosynthesis from taking place until the 
temperature drops below 100 keV.) At that point the neutron/proton ratio is approximately 
1/7. Of all the light nuclei, it is energetically favorable for the nucleons to reside in ^He, 
and indeed that is what most of the free neutrons are converted into; for every two neutrons 
and fourteen protons, we end up with one helium nucleus and twelve protons. Thus, about 
25% of the baryons by mass are converted to helium. In addition, there are trace amounts of 
deuterium (approximately 10~^ deuterons per proton), ^He (also ~ 10~^), and ^Li (~ 10"^'^). 

Of course these numbers are predictions, which are borne out by observations of the 
primordial abundances of light elements. (Heavier elements are not synthesized in the Big 
Bang, but require supernova explosions in the later universe.) We have glossed over numerous 
crucial details, especially those which explain how the different abundances depend on the 
cosmological parameters. For example, imagine that we deviate from the Standard Model by 
introducing more than three light neutrino species. This would increase the radiation energy 
density at a fixed temperature, which in turn decreases the timescales associated with a given 
temperature (since t ~ oc Pr^^^). Nucleosynthesis would therefore happen somewhat 
earlier, resulting in a higher abundance of neutrons, and hence in a larger abundance of ^He. 
Observations of the primordial helium abundance, which arc consistent with the Standard 
Model prediction, provided the first evidence that the number of light neutrinos is actually 
three. 

The most amazing fact about nucleosynthesis is that, given that the universe is radiation 
dominated during the relevant epoch (and the physics of general relativity and the Stan- 
dard Model), the relative abundances of the light elements depend essentially on just one 
parameter, the baryon to entropy ratio 

_ ub n-ni 



where ub — — ni \s the difference between the number of baryons and antibaryons per 
unit volume. The range of 77 consistent with the deuterium and ^He primordial abundances 
is 

2.6 X 10"^° < 77 < 6.2 X 10"^° . (127) 
Very recently this number has been independently determined to be 

= 6.1 X 10-° t^ll^^Z (128) 

from precise measurements of the relative heights of the first two microwave background 
(CMB) acoustic peaks by the WMAP satellite. This is illustrated in figure (4.9) and we will 
have a lot more to say about this quantity later when we discuss baryogenesis. 
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Figure 4.9: Abundances of light elements produced by BBN, as a function of the baryon den- 
sity. The vertical strip indicates the concordance region favored by observations of primordial 
abundances. From [47]. 

4.6 Finite Temperature Phase Transitions 

We have hinted at the need to go beyond perfect fluid sources for the Einstein equations if 
we are to unravel some of the mysteries left by the standard cosmology. Given our modern 
understanding of particle physics, it is natural to consider using field theory to model matter 
at early times in the universe. The effect of cosmic expansion and the associated thermody- 
namics yield some fascinating phenomena when combined with such a field theory approach. 
One significant example is provided by finite temperature phase transitions. 

Rather than performing a detailed calculation in finite-temperature field theory, we will 
illustrate this with a rough argument. Consider a theory of a single real scalar field at zero 
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temperature, interacting with a second real scalar field x- The Lagrangian density is 



c = -\d,<pd^<p - \d,xd^x - v{<p,T = 0) - I^V , (129) 

with potential 

y((/.,T = 0) = -^(/.^ + ^</.\ (130) 

where is a parameter with dimensions of mass and A and g are dimensionless coupling 
constants. 

Now, consider what the effective theory for looks like if we assume that x is in thermal 
equilibrium. In this case, we may replace x by the temperature T to obtain a Lagrangian 
for (j) only, with finite temperature effective potential given by 

V{<P,T) = \{gT^-^^)<P^ + \ct>^ . (131) 

This simple example demonstrates a very significant result. At zero temperature, all your 
particle physics intuition tells you correctly that the theory is spontaneously broken, in this 
case the global Z2 symmetry of the Lagrangian is broken by the ground state value (0) 7^ 
(we shall see more of this soon). However, at temperatures above the critical temperature 
Tc given by 

7; = -^, (132) 

the full Z2 symmetry of the Lagrangian is respected by the finite temperature ground state, 
which is now at (0) = 0. This behavior is know as finite temperature symmetry restoration, 
and its inverse, occurring as the universe cools, is called finite temperature spontaneous 
symmetry breaking. 

Zero temperature unification and symmetry breaking are fundamental features of modern 

particle physics models. For particle physicists this amounts to the statement that the chosen 
vacuum state of a gauge field theory is one which does not respect the underlying symmetry 
group of the Lagrangian. 

The melding of ideas from particle physics and cosmology described above leads us to 
speculate that matter in the early universe was described by a unified gauge field theory based 
on a simple continuous Lie group, G. As a result of the extreme temperatures of the early 
universe, the vacuum state of the theory respected the full symmetry of the Lagrangian. As 
the universe cooled, we hypothesize that the gauge theory underwent a series of spontaneous 
symmetry breakings (SSB) until matter was finally described by the unbroken gauge groups 
of QCD and QED. Schematically the symmetry breaking can be represented by 

G^H^ . 5t/(3)e X SU{2)l X U{1)y ^ -St/(3)e x U{l),m ■ (133) 

The group G is known as the grand unified gauge group and the initial breaking G ^ H 
and is expected to take place at lO^^GeV, as we discuss below. 
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Let us briefly set up the mathematical description of SSB which will be useful later when 
we discuss topological properties of the theory. 

Consider a gauge field theory described by a continuous group G. Denote the vacuum 
state of the theory by |0) . Then, given g{x) e G, the state g{x)\0) is also a vacuum. Suppose 
all vacuum states are of the form (7|0) but that the state |0) is invariant only under a subgroup 
H G G. Then the symmetry of the theory is said to have been spontaneously broken from G 
to H. Let us define two vacuum states |0)^ and \0)b to represent the same state under the 
broken group if g\0)A = |0)b for some g E G. We write |0)^ ~ \0)b and the distinguishable 
vacua of the theory are equivalence classes under ~ and are the cosets of if in G. Then the 
vacuum manifold, the space of all accessible vacua of the theory, is the coset space 

M = G/H . (134) 

Given that we believe the universe evolved through a sequence of such symmetry breakings, 
there are many questions we can ask about the cosmological implications of the scheme and 
many ways in which we can constrain and utilize the possible breaking patterns. 

4.7 Topological Defects 

In a quantum field theory, small oscillations around the vacuum appear as particles. If the 
space of possible vacuum states is topologically nontrivial, however, there arises the possibil- 
ity of another kind of solitonic object: a topological defect. Consider a field theory described 
by a continuous symmetry group G which is spontaneously broken to a subgroup H d G. 
Recall that the space of all accessible vacua of the theory, the vacuum manifold, is defined 
to be the space of cosets oi H m G\ M. = G/H. Whether the theory admits topological 
defects depends on whether the vacuum manifold has nontrivial homotopy groups. A homo- 
topy group consists of equivalence classes of maps of spheres (with fixed base point) into the 
manifold, where two maps are equivalent if they can be smoothly deformed into each other. 
The homotopy groups defined in terms of n-spheres are denoted 7r„. 

In general, a field theory with vacuum manifold Jv[ possesses a topological defect of some 
type if 

T^i{M) 7^ 1 , (135) 

for some i = 0, 1, . . .. In particular, we can have a set of defects, as listed in table 3. In 
order to get an intuitive picture of the meaning of these topological criteria let us consider 
the first two. 

If 'Kq{M) 7^ 1, the manifold M. is disconnected. (A zero-sphere is the set of two points a 
fixed distance from the origin in R^. The set of topologically equivalent maps with fixed base 
point from such such a sphere into a manifold is simply the set of disconnected pieces into 
which the manifold falls.) Let's assume for simplicity that the vacuum manifold consists 
of just two disconnected components A^i and M.2 and restrict ourselves to one spatial 
dimension. Then, if we apply boundary conditions that the vacuum at — oo lies in M.i and 
that at -|-oo lies in A^2, by continuity there must be a point somewhere where the order 
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Homotopy Constraint 


Topological Defect 


MM) ^ 1 


Domain Wall 


MM) ^ 1 


Cosmic String 


MM) ^ 1 


Monopole 


MM) ^ 1 


Texture 



Table 3: Topological defects, as described by homotopy groups of the vacuum manifolds of 
a theory with spontaneously broken symmetry. 

parameter does not lie in the vacuum manifold as the field interpolates between the two 
vacua. The region in which the field is out of the vacuum is known as a domain wall. 

Similarly, suppose 7ri(A^) ^ 1. This implies that the vacuum manifold is not simply- 
connected: there are non-contractible loops in the manifold. The order parameter for such 
a situation may be considered complex and if, when traveling around a closed curve in 
space, the phase of the order parameter changes by a non-zero multiple of 2tt, by continuity 
there is a point within the curve where the field is out of the vacuum manifold. Continuity 
in the direction perpendicular to the plane of the curve implies that there exists a line of 
such points. This is an example of a cosmic string or line defect. Since we shall mostly 
concentrate on these defects let us now give a detailed analysis of the simplest example, the 
Nielsen- Olesen vortex in the Ahelian Higgs model. 

Consider a complex scalar field theory based on the Abelian gauge group U{1). The 
Lagrangian density for this model is 

C = ~{D,^)*D^^ - \f,^F^^ - V{^) , (136) 

where ^p{x) is the complex scalar field, the covariant derivative, D^, and field strength tensor 
F^y are defined in terms of the Abelian gauge field as 

L>^(^ = {d^ + ieA^,)^p , (137) 

F^u = d^A, - d,A^ , (138) 
and the symmetry breaking (or "Mexican hat" ) potential is 

y(9.) = ^(^V-r/Y. (139) 

Here e is the gauge coupling constant and is a parameter that represents the scale of the 
symmetry breaking. In preparation for our discussion of the more general cosmological situ- 
ation, let us note that at high temperatures the Lagrangian contains temperature-dependent 
corrections given by 

v{^p) ^ v{^) + crV' + • • • , (140) 
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where C is a constant. For T > Tc, where is defined by 

" 2C 



= ^rf , (141) 



we see that 



V{^) = ^((^V) + ^r}'{-^l]^*^+ '-v' (142) 




is minimized by {ip) — 0. However, for T < Tc, minimizing the above expression with respect 
to (p yields a minimum at 

-V'{l-^) . (143) 

Thus, for T > Tc the full symmetry group G = U{1) is restored and the vacuum expectation 
value of the 99-field is zero. As the system cools, there is a phase transition at the critical 
temperature Tc and the vacuum symmetry is spontaneously broken, in this case entirely: 

G = U{1) — >1 = H . (144) 

The vacuum manifold is 

M = G/H = U{l)/l = U(l) = {(p : (/? = rje''' , < a < 27r} . (145) 

The group U{1) is topologically a circle, so 

M^SK (146) 

The set of topologically equivalent ways to map one circle to another circle is given by the 
integers. (We can wrap it any number of times in either sense.) The first homotopy group 
is therefore the integers. Calculating homotopy groups in general is difficult, but for all 
of the other groups are trivial. We therefore have 

MS') = 1 (147) 

7ri(5^) = Z (148) 

MS') = 1 (149) 

MS') = 1- (150) 

The Abelian Higgs model therefore allows for cosmic strings. 

To see how strings form in this model, consider what happens as T decreases through Tc. 
The symmetry breaks and the field acquires a VEV given by 

{if) = 776- , (151) 

where a may be chosen differently in different regions of space, as implied by causality (we 
shall discuss this shortly). The requirement that (ip) be single valued implies that around 
any closed curve in space the change Aa in a must satisfy 

Aa = 27rn , nE Z . (152) 
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If for a given loop we have n 7^ 0, then we can see that any 2-surface bounded by the loop 
must contain a singular point, for if not then we can continuously contract the loop to a 
point, implying that n — which is a contradiction. At this singular point the phase, a, 
is undefined and ((/?) = 0. Further, {(f) must be zero all along an infinite or closed curve, 
since otherwise wc can contract our loop without encountering a singularity. We identify 
this infinite or closed curve of false vacuum points as the core of our string. We shall restrict 
our attention to the case where n = 1, since this is the most likely configuration. 

The first discussion of string solutions to the Abelian Higgs model is due to Nielsen 
and Olesen. Assume the string is straight and that the core is aligned with the 2;- axis. In 
cylindrical polar coordinates, (r, 9), let us make the ansatz 



= ^Yir)d^e. (153) 
The Lagrangian then yields the simplified equations of motion 

-X"-^ + ^ + ^X{X--l) = 0, 
r 2 

+ ^ + ^ 0, (154) 

r A 

where a prime denotes differentiation with respect to r. It is not possible to solve these 
equations analytically, but asymptotically we have 

X r^r , y = 1 + 0{r^) ; as r ^ , 
/ -aVv\ 

Xr^l-Ol , y ~ 0(^e-^ ' ; as r ^ 00 . (155) 

This solution corresponds to a string centered on the ;2-axis, with a central magnetic core of 
width ~ {erj)~^ carrying a total magnetic flux 

$ = i A^dx*" = . (156) 

Jr=oo 6 

The core region over which the Higgs fields are appreciably non-zero has width ~ (A^^^?])"^. 
Note that these properties depend crucially on the fact that we started with a gauged sym- 
metry. A spontaneously broken global symmetry with nontrivial Tii{M) would still produce 
cosmic strings, but they would be much less localized, since there would be no gauge fields 
to cancel the scalar gradients at large distances. 

Strings are characterized by their tension^ which is the energy per unit length. For the 
Nielsen-Olesen solution just discussed, the tension is approximately 

^lr^rf . (157) 
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Thus, the energy of the string is set by the expectation value of the order parameter respon- 
sible for the symmetry breaking. This behavior is similarly characteristic of other kinds of 
topological defects. Again, global defects are quite different; the tension of a global string 
actually diverges, due to the slow fall-off of the energy density as we move away from the 
string core. It is often convenient to parameterize the tension by the dimensionless quantity 
Gfi, so that a Planck-scale string would have Gfi ~ 1. 

From the above discussion we can see that cosmology provides us with a unique op- 
portunity to explore the rich and complex structure of particle physics theories. Although 
the topological solutions discussed above exist in the theory at zero temperature, there is 
no mechanism within the theory to produce these objects. The topological structures con- 
tribute a set of zero measure in the phase space of possible solutions to the theory and hence 
the probability of production in particle processes is exponentially suppressed. What the 
cosmological evolution of the vacuum supplies is a concrete causal mechanism for producing 
these long lived exotic solutions. 

At this point it is appropriate to discuss how many of these defects we expect to be 
produced at a cosmological phase transition. In the cosmological context, the mechanism 
for the production of defects is known as the Kibble mechanism. The guiding principle here 
is causality. As the phase transition takes place, the maximum causal distance imposed on 
the theory by cosmology is simply the Hubble distance - the distance which light can have 
traveled since the big bang. 

As we remarked above, as the temperature of the universe falls well below the critical 
temperature of the phase transition, Tc, the expectation value of the order parameter takes 
on a definite value (~ rye*" in the Abelian Higgs model) in each region of space. However, 
at temperatures around the critical temperature we expect that thermal fluctuations in (0) 
will be large so that as the universe cools it will split into domains with different values of 
a in different domains. This is the crucial role played by the cosmological evolution. 

The Kibble mechanism provides us with an order of magnitude upper bound for the 
size of such a region as the causal horizon size at the time of the phase transition. The 
boundaries between the domains will be regions where the phase of (0) changes smoothly. 
If the phase changes by 27m for some n ^ when traversing a loop in space, then any 
surface bounded by that loop intersects a cosmic string. These strings must be horizon-sized 
or closed loops. Numerical simulations of cosmic string formation indicate that the initial 
distribution of strings consists of 80% horizon-sized and 20% loops by mass. If we consider 
the temperature at which there is insufficient thermal energy to excite a correlation volume 
back into the unbroken state, the Ginsburg temperature, Tq, then, given the assumption of 
thermal equilibrium above the phase transition, a much improved estimate for the initial 
separation of the defects can be derived and is given by 

ate) ~ A-^-^ , 

where A is the self coupling of the order parameter. This separation is microscopic. 

Interesting bounds on the tension of cosmic strings produced by the Kibble mechanism 
arise from two sources: perturbations of the CMB, and gravitational waves. Both arise 
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because the motions of heavy strings moving at relativistic velocities lead to time-dependent 
gravitational fields. The actual values of the bounds are controversial, simply because it is 
difficult to accurately model the nonlinear evolution of a string network, and the results can 
be sensitive to what assumptions are made. Nevertheless, the CMB bounds amount roughly 
to [118, 119, 120, 121] 

Gfx < 10"^ . (158) 

This corresponds roughly to strings at the GUT scale, 10^^ GeV, which is certainly an 
interesting value. Bounds from gravitational waves come from two different techniques: 
direct observation, and indirect measurement through accurate pulsar timings. Currently, 
pulsar timing measurements are more constraining, and give a bound similar to that from 
the CMB [122] . Unlike the CMB measurements, however, gravitational wave observatories 
will become dramatically better in the near future, through operations of ground-based 
observatories such as LIGO as well as satellites such as LISA. These experiments should be 
able to improve the bounds on G/i by several orders of magnitude [123]. 

Now that we have established the criteria necessary for the production of topological 
defects in spontaneously broken theories, let us apply these conditions to the best understood 
physical example, the electroweak phase transition. As we remarked earlier, the GWS theory 
is based on the gauge group SU{2)l x f/(l)y. At the phase transition this breaks to pure 
electromagnetism, U{l)em- Thus, the vacuum manifold is the space of cosets 

Mew = [SU(2) x U(l)]/U(l) , (159) 

This looks complicated but in fact this space is topologically equivalent to the three-sphere, 
S^. The homotopy groups of the three-sphere are 

MS') = 1 (160) 

TTiiS^) = 1 (161) 

MS') = 1 (162) 

MS') = (163) 

(Don't be fooled into thinking that all homotopy groups of spheres vanish except in the 
dimensionality of the sphere itself; for example, 7r3(S'^) = Z.) 

Thus, the electroweak model does not lead to walls, strings, or monopoles. It does lead 
to what we called "texture," which deserves further comment. In a theory where 7r3(A1) is 
nontrivial but the other groups vanish, we can always map three-dimensional space smoothly 
into the vacuum manifold; there will not be a defect where the field chmbs out of M.. How- 
ever, if we consider field configurations which approach a unique value at spatial infinity, 
they will fall into homotopy classes characterized by elements of nzi^hA)] configurations with 
nonzero winding will be textures. If the symmetry is global, such configurations will nec- 
essarily contain gradient energies from the scalar fields. The energy perturbations caused 
by global textures were, along with cosmic strings, formerly popular as a possible origin of 



52 



structure formation in the universe [124, 125]; the predictions of these theories are incon- 
sistent with the sharp acoustic peaks observed in the CMB, so such models are no longer 
considered viable. 

In the standard model, however, the broken symmetry is gauged. In this case there is 
no need for gradient energies, since the gauge field can always be chosen to cancel them; 
equivalcntly, "texture" configurations can always be brought to the vacuum by a gauge 
transformation. However, transitions from one texture configuration to one with a different 
winding number are gauge invariant. These transitions will play a role in electroweak baryon 
number violation, discussed in the next section. 



4.8 Baryogenesis 

The symmetry between particles and antiparticles [126, 127], firmly established in coUider 
physics, naturally leads to the question of why the observed universe is composed almost 
entirely of matter with little or no primordial antimatter. 

Outside of particle accelerators, antimatter can be seen in cosmic rays in the form of a 
few antiprotons, present at a level of around 10"*^ in comparison with the number of protons 
(for example see [128]). However, this proportion is consistent with secondary antiproton 
production through accelerator-like processes, p -|- p — > 3p + p, as the cosmic rays stream 
towards us. Thus there is no evidence for primordial antimatter in our galaxy. Also, if matter 
and antimatter galaxies were to coexist in clusters of galaxies, then we would expect there 
to be a detectable background of 7-radiation from nucleon-antinucleon annihilations within 
the clusters. This background is not observed and so we conclude that there is negligible 
antimatter on the scale of clusters (For a review of the evidence for a baryon asymmetry see 
[129].) 

More generally, if large domains of matter and antimatter exist, then annihilations would 
take place at the interfaces between them. If the typical size of such a domain was small 
enough, then the energy released by these annihilations would result in a diffuse 7-ray back- 
ground and a distortion of the cosmic microwave radiation, neither of which is observed. 

While the above considerations put an experimental upper bound on the amount of an- 
timatter in the universe, strict quantitative estimates of the relative abundances of baryonic 
matter and antimatter may also be obtained from the standard cosmology. The baryon num- 
ber density does not remain constant during the evolution of the universe, instead scaling 
like a~^, where a is the cosmological scale factor [130]. It is therefore convenient to define 
the baryon asymmetry of the universe in terms of the quantity 

V^'f, (164) 

defined earlier. Recall that the range of 7] consistent with the deuterium and ^He primordial 
abundances is 

2.6 X 10-^° < 77 < 6.2 X 10-^° . (165) 

Thus the natural question arises; as the universe cooled from early times to today, what 
processes, both particle physics and cosmological, were responsible for the generation of 
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this very specific baryon asymmetry? (For reviews of meclianisms to generate tlie baryon 
asymmetry, see 

As pointed out by Sakfiarov [131], a small baryon asymmetry r) may have been produced 
in the early universe if three necessary conditions are satisfied 

• baryon number {B) violation, 

• violation of C (charge conjugation symmetry) and CP (the composition of parity and 



• departure from thermal equilibrium. 

The first condition should be clear since, starting from a baryon symmetric universe with 
?7 = 0, baryon number violation must take place in order to evolve into a universe in which 77 
does not vanish. The second Sakharov criterion is required because, if C and CP are exact 
symmetries, one can prove that the total rate for any process which produces an excess 
of baryons is equal to the rate of the complementary process which produces an excess of 
antibaryons and so no net baryon number can be created. That is to say that the thermal 
average of the baryon number operator B, which is odd under both C and CP, is zero unless 
those discrete symmetries are violated. CP violation is present cither if there arc complex 
phases in the Lagrangian which cannot be reabsorbed by field redefinitions (explicit breaking) 
or if some Higgs scalar field acquires a VEV which is not real (spontaneous breaking). We 
will discuss this in detail shortly. 

Finally, to explain the third criterion, one can calculate the equilibrium average of B at 
a temperature T — 1//3: 



where we have used that the Hamiltonian H commutes with CPT. Thus {B)t = in 
equilibrium and there is no generation of net baryon number. 

Of the three Sakharov conditions, baryon number violation and C and CP violation may 
be investigated only within a given particle physics model, while the third condition - the 
departure from thermal equilibrium - may be discussed in a more general way, as we shall 
see (for baryogenesis reviews see [132, 133, 9, 10, 134, 135, 136].) Let us discuss the Sakharov 
criteria in more detail. 

4.9 Baryon Number Violation 

4.9.1 violation in Grand Unified Theories 

As discussed earher. Grand Unified Theories (GUTs) [137] describe the fundamental inter- 
actions by means of a unique gauge group G which contains the Standard Model (SM) gauge 
group SU (3)c ® SU {2)l®U (l)y. The fundamental idea of GUTs is that at energies higher 



C), 



{B)t 



Tr {e-^^^'B) = Tr [{C PT){C PT)-\-^^ B)] 
Tr {e-^"{CPT)-^B{CPT)\ = -Tr (e^^^E) , 



(166) 
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than a certain energy threshold Mqut the group symmetry is G and that, at lower ener- 
gies, the symmetry is broken down to the SM gauge symmetry, possibly through a chain of 
symmetry breakings. The main motivation for this scenario is that, at least in supcrsym- 
metric models, the (running) gauge couphngs of the SM unify [138, 139, 140] at the scale 
Mgut ~ 2 X 10^^ GeV, hinting at the presence of a GUT involving a higher symmetry with 
a single gauge coupling. 

Baryon number violation seems very natural in GUTs. Indeed, a general property of 
these theories is that the same representation of G may contain both quarks and leptons, 
and therefore it is possible for scalar and gauge bosons to mediate gauge interactions among 
fermions having different baryon number. 

4.9.2 5-violation in the Electroweak theory. 

It is well-known that the most general renormalizable Lagrangian invariant under the SM 
gauge group and containing only color singlet Higgs fields is automatically invariant under 
global abelian symmetries which may be identified with the baryonic and leptonic symme- 
tries. These, therefore, are accidental symmetries and as a result it is not possible to violate 
B and L at tree-level or at any order of perturbation theory. Nevertheless, in many cases 
the perturbative expansion does not describe all the dynamics of the theory and, indeed, 
in 1976 't Hooft [141] realized that nonperturbative effects (instantons) may give rise to 
processes which violate the combination B -\- L, but not the orthogonal combination B — L. 
The probability of these processes occurring today is exponentially suppressed and probably 
irrelevant. However, in more extreme situations - like the primordial universe at very high 
temperatures [142, 143, 144, 145] - baryon and lepton number violating processes may be 
fast enough to play a significant role in baryogenesis. Let us have a closer look. 

At the quantum level, the baryon and the lepton symmetries are anomalous [146, 147], 
so that their respective Noether currents and jj^ are no longer conserved, but satisfy 

d,fB = d,fL = nf y^,W;^W'^^^ - ^.F.^F^"^ , (167) 

where g and g' are the gauge couplings of SU{2)l and f/(l)y, respectively, nj is the number 
of families and W^" = {l/2)e^^''"-^Wap is the dual of the SU{2)l field strength tensor, with 
an analogous expression holding for F. To understand how the anomaly is closely related 
to the vacuum structure of the theory, we may compute the change in baryon number from 
time i = to some arbitrary final time t = tf. For transitions between vacua, the average 
values of the field strengths are zero at the beginning and the end of the evolution. The 
change in baryon number may be written as 

AB = ANcs = nf[Ncs{tf) - NcsiO)] . (168) 
where the Chern-Simons number is defined to be 

Ncs{t) = d'xe'^'^ Tr (^A^djA^ + '^igAAjA^^ , (169) 
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where Ai is the 811(2)^ gauge field. Although the Chern-Simons number is not gauge in- 
variant, the change ANqs is. Thus, changes in Chern-Simons number result in changes in 
baryon number which are integral multiples of the number of families n/ (with n/ = 3 in the 
real world). Gauge transformations U{x) which connects two degenerate vacua of the gauge 
theory may change the Chern-Simons number by an integer n, the winding number. 

If the system is able to perform a transition from the vacuum Q^^l ^he closest one 
the Chern-Simons number is changed by unity and Ai? = AL = rif. Each transition 
creates 9 left-handed quarks (3 color states for each generation) and 3 left-handed leptons 
(one per generation). 

However, adjacent vacua of the electroweak theory are separated by a ridge of configu- 
rations with energies larger than that of the vacuum. The lowest energy point on this ridge 
is a saddle point solution to the equations of motion with a single negative eigenvalue, and 
is referred to as the sphaleron [143, 144]. The probability of baryon number nonconserving 
processes at zero temperature has been computed by 't Hooft [141] and is highly suppressed. 

The thermal rate of baryon number violation in the broken phase is 



where /x is a dimensionless constant. Although the Boltzmann suppression in (170) appears 
large, it is to be expected that, when the electroweak symmetry becomes restored at a 
temperature of around 100 GeV, there will no longer be an exponential suppression factor. 
A simple estimate is that the rate per unit volume of sphaleron events is 



with K another dimensionless constant. The rate of sphaleron processes can be related to 
the diffusion constant for Chern-Simons number by a fiuctuation-dissipation theorem [148] 
(for a good description of this see [134]). 

4.9.3 CP violation 

CP violation in GUTs arises in loop-diagram corrections to baryon number violating bosonic 
decays. Since it is necessary that the particles in the loop also undergo S-violating decays, 
the relevant particles are the X, Y, and bosons in the case of SU{5). 

In the electroweak theory things are somewhat different. Since only the left-handed 
fermions are SU{2)l gauge coupled, C is maximally broken in the SM. Moreover, CP is 
known not to be an exact symmetry of the weak interactions. This is seen experimentally 
in the neutral kaon system through Kq, Kq mixing. Thus, CP violation is a natural feature 
of the standard electroweak model. 

While this is encouraging for baryogenesis, it turns out that this particular source of CP 
violation is not strong enough. The relevant effects are parameterized by a dimensionless 
constant which is no larger than lO"^'^. This appears to be much too small to account for the 
observed BAU and, thus far, attempts to utihze this source of CP violation for electroweak 




(170) 



(171) 
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baryogenesis have been unsuccessful. In light of this, it is usual to extend the SM in some 
fashion that increases the amount oi CP violation in the theory while not leading to results 
that conflict with current experimental data. One concrete example of a well-motivated 
extension in the minimal supersymmetric standard model (MSSM). 

4.9.4 Depcirture from Thermal Equilibrium 

In some scenarios, such as GUT baryogenesis, the third Sakharov condition is satisfied due 
to the presence of superheavy decaying particles in a rapidly expanding universe. These 
generically fall under the name of out-of-equilibrium decay mechanisms. 

The underlying idea is fairly simple. If the decay rate Tx of the superheavy particles X 
at the time they become nonrelativistic {i.e. at the temperature T ~ Mx) is much smaller 
than the expansion rate of the universe, then the X particles cannot decay on the time scale 
of the expansion and so they remain as abundant as photons for T ^ Mx- In other words, 
at some temperature T > Mx, the superheavy particles X are so weakly interacting that 
they cannot catch up with the expansion of the universe and they decouple from the thermal 
bath while they are still relativistic, so that nx ^ at the time of decoupling. 

Therefore, at temperature T ~ Mx, they populate the universe with an abundance which 
is much larger than the equilibrium one. This overabundance is precisely the departure from 
thermal equilibrium needed to produce a final nonvanishing baryon asymmetry when the 
heavy states X undergo B and CP violating decays. 

The out-of-equihbrium condition requires very heavy states: Mx (10^^ — 10^^) GeV 
and Mx ^ (10^° — 10^^) GeV, for gauge and scalar bosons, respectively [149], if these heavy 
particles decay through rcnormalizablc operators. 

A different implementation can be found in the electroweak theory. At temperatures 
around the electroweak scale, the expansion rate of the universe in thermal units is small 
compared to the rate of baryon number violating processes. This means that the equilibrium 
description of particle phenomena is extremely accurate at electroweak temperatures. Thus, 
baryogenesis cannot occur at such low scales without the aid of phase transitions and the 
question of the order of the electroweak phase transition becomes central. 

If the EWPT is second order or a continuous crossover, the associated departure from 
equilibrium is insufficient to lead to relevant baryon number production [145]. This means 
that for EWBG to succeed, we either need the EWPT to be strongly first order or other 
methods of destroying thermal equilibriTim to be present at the phase transition. For a first 
order transition there is an extremum at = which becomes separated from a second local 
minimum by an energy barrier. At the critical temperature T = Tc both phases are equally 
favored energetically and at later times the minimum at ip becomes the global minimum 
of the theory. 

The dynamics of the phase transition in this situation is crucial to most scenarios of 
electroweak baryogenesis. The essential picture is that around Tc quantum tunneling occurs 
and nucleation of bubbles of the true vacuum in the sea of false begins. Initially these bubbles 
are not large enough for their volume energy to overcome the competing surface tension and 
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they shrink and disappear. However, at a particular temperature below T^, bubbles just large 
enough to grow nucleate. These are termed critical bubbles, and they expand, eventually 
filling all of space and completing the transition. 

As the bubble walls pass each point in space, the order parameter changes rapidly, as do 
the other fields, and this leads to a significant departure from thermal equilibrium. Thus, if 
the phase transition is strongly enough first order it is possible to satisfy the third Sakharov 
criterion in this way. 

A further natural way to depart from equilibrium is provided by the dynamics of topo- 
logical defects [150, 151, 152, 153]. If, for example, cosmic strings are produced at the GUT 
phase transition, then the decays of loops of string act as an additional source of superheavy 
bosons which undergo baryon number violating decays. 

When defects are produced at the TeV scale, a detailed analysis of the dynamics of a 
network of these objects shows that a significant baryon to entropy ratio can be generated if 
the electroweak symmetry is restored around such a higher scale ordinary defect. Although 
S-violation can be inefficient along nonsuperconducting strings [154], there remain viable 
scenarios involving other ordinary defects, superconducting strings or defects carrying baryon 
number. 

4.9.5 Baryogenesis via leptogenesis 

Since the linear combination B — L is left unchanged by sphaleron transitions, the baryon 
asymmetry may be generated from a lepton asymmetry [155, 156, 157] (see also [158, 159].) 
Indeed, sphaleron transition will reprocess any lepton asymmetry and convert (a fraction 
of) it into baryon number. This is because B + L must be vanishing and the final baryon 
asymmetry results to be 5 ~ —L. 

In the SM as well as in its unified extension based on the group SU (5) , 5 — L is conserved 
and no asymmetry in B — L can be generated. However, adding right-handed Majorana neu- 
trinos to the SM breaks B — L and the primordial lepton asymmetry may be generated by the 
out-of-equilibrium decay of heavy right-handed Majorana neutrinos (in the supersym- 
metric version, heavy scalar neutrino decays are also relevant for leptogenesis). This simple 
extension of the SM can be embedded into GUTs with gauge groups containing SO (10). 
Heavy right-handed Majorana neutrinos can also explain the smallness of the light neutrino 
masses via the see-saw mechanism [157]. 

4.9.6 Affleck-Dine Beiryogenesis 

Finally in this section, we mention briefiy a mechanism introduced by Affleck and Dine [160] 
involving the cosmological evolution of scalar fields carrying baryonic charge. 

Consider a colorless, electrically neutral combination of quark and lepton fields. In a su- 
pcrsymmetric theory this object has a scalar superpartner, %, composed of the corresponding 
squark q and slepton I fields. 

An important feature of supersymmetric field theories is the existence of "fiat directions" 
in field space, on which the scalar potential vanishes. Consider the case where some com- 
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ponent of the field x li^s along a flat direction. By this we mean that there exist directions 
in the superpotential along which the relevant components of x can be considered as a free 
massless field. At the level of renormalizable terms, fiat directions are generic, but super- 
symmetry breaking and nonrenormalizable operators lift the fiat directions and sets the scale 
for their potential. 

During inflation it is likely that the x field is displaced from the position (x) = 0, 
establishing the initial conditions for the subsequent evolution of the field. An important 
role is played at this stage by baryon number violating operators in the potential V^(x), 
which determine the initial phase of the field. When the Hubble rate becomes of the order 
of the curvature of the potential ~ "^3/2, the condensate starts oscillating around its present 
minimum. At this time, 5-violating terms in the potential are of comparable importance 
to the mass term, thereby imparting a substantial baryon number to the condensate. After 
this time, the baryon number violating operators are negligible so that, when the baryonic 
charge of x is transferred to fermions through decays, the net baryon number of the universe 
is preserved by the subsequent cosmological evolution. 

The most recent implementations of the Afiicck-Dine scenario have been in the context 
of the minimal supersymmetric standard model [161, 162] in which, because there are large 
numbers of fields, flat directions occur because of accidental degeneracies in field space. 



5 Inflation 

In our previous lectures we have described what is known as the standard cosmology. This 
framework is a towering achievement, describing to great accuracy the physical processes 
leading to the present day universe. However, there remain outstanding issues in cosmology. 
Many of these come under the heading of initial condition problems and require a more 
complete description of the sources of energy density in the universe. The most severe of 
these problems eventually led to a radical new picture of the physics of the early universe - 
cosmological inflation [163, 164, 165], which is the subject of this lecture. 

We will begin by describing some of the problems of the standard cosmology. 



5.1 The Flatness Problem 

The Friedmann equation may be written as 



^-l-THTi' (172) 



k_ 

where for brevity we are now writing Q instead of fitotai- Differentiating this with respect to 
the scale factor, this imphes 

— - l + 3w ^ ^. 173 

da a 

This equation is easily solved, but its most general properties are all that we shall need and 
they are qualitatively different depending on the sign of 1 -|-3u'. There are three fixed points 
of this differential equation, as given in table 4. 
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Fixed Point 


1 + 3w > 


1 + 3w < 


Q = 


attractor 


repeller 


fi = 1 


repeller 


attractor 


Q, — oo 


attractor 


repeller 



Table 4: Behavior of the density parameter near fixed points. 
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Figure 5.10: Past light cones in a universe expanding from a Big Bang singularity, illustrating 
particle horizons in cosmology. Points at recombination, observed today as parts of the 
cosmic microwave background on opposite sides of the sky, have non-overlapping past light 
cones (in conventional cosmology); no causal signal could have infiuenced them to have the 
same temperature. 

Obscrvationally we know that ~ 1 today - i.e., we are very close to the repeller of this 
differential equation for a universe dominated by ordinary matter and radiation [w > — 1/3) . 
Even if we only took account of the luminous matter in the universe, we would clearly live in 
a universe that was far from the attractor points of the equation. It is already quite puzzling 
that the universe has not reached one of its attractor points, given that the universe has 
evolved for such a long time. However, we may be more quantitative about this. If the only 
matter in the universe is radiation and dust, then in order to have fl in the range observed 
today requires (conservatively) 

< 1 - Q < 10-60 . (174) 

This remarkable degree of fine tuning is the fiatness problem. Within the context of the 
standard cosmology there is no known explanation of this fine-tuning. 

5.2 The Horizon Problem 

The horizon problem stems from the existence of particle horizons in FRW cosmologies, 
as discussed in the first lecture. Horizons exist because there is only a finite amount of 
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time since the Big Bang singularity, and thus only a finite distance that photons can travel 
within the age of the universe. Consider a photon moving along a radial trajectory in a flat 
universe (the generalization to non-flat universes is straightforward). In a flat universe, we 
can normalize the scale factor to 

ao = 1 (175) 
without loss of generality. A radial null path obeys 

Q = ds^ = -df + a^dr^ , (176) 

so the comoving (coordinate) distance traveled by such a photon between times ti and t2 is 

Ar -r4,. (177) 
iti a(t) ^ ' 

To get the physical distance as it would be measured by an observer at any time t, simply 
multiply by a(t). For simplicity let's imagine we are in a matter-dominated universe, for 
which 

The Hubble parameter is therefore given by 



- ■ (178) 



H = Hr' 

3 

= a-3/2i/o . (179) 

Then the photon travels a comoving distance 

Ar = 2i/o-i ( - ^) . (180) 

The comoving horizon size when a = a* is the distance a photon travels since the Big Bang, 

rhor(a*) = 2i/o-iy^ . (181) 

The physical horizon size, as measured on the spatial hypersurface at a*, is therefore simply 

4or(a*) = a*nior(a*) = 2-^^"^ . (182) 

Indeed, for any nearly-flat universe containing a mixture of matter and radiation, at any one 
epoch we will have 

4or(a*) ~ i^-' , (183) 

where is the Hubble distance at that particular epoch. This approximate equality 

leads to a strong temptation to use the terms "horizon distance" and "Hubble distance" 
interchangeably; this temptation should be resisted, since inflation can render the former 
much larger than the latter, as we will soon demonstrate. 
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The horizon problem is simply the fact that the CMB is isotropic to a high degree of 
precision, even though widely separated points on the last scattering surface are completely 
outside each others' horizons. When we look at the CMB we were observing the universe 
at a scale factor ocmb ~ 1/1200; meanwhile, the comoving distance between a point on the 
CMB and an observer on Earth is 

Ar = 2i7o"' (1 - y^a^) 

Ri 2Hq^ . (184) 

However, the comoving horizon distance for such a point is 

^hor(aCMB) = 277,7 ^V'^CMB 

6 X lO-^i^o"^ ■ (185) 

Hence, if we observe two widely-separated parts of the CMB, they will have non-overlapping 
horizons; distinct patches of the CMB sky were causally disconnected at recombination. 

Nevertheless, they are observed to be at the same temperature to high precision. The 
question then is, how did they know ahead of time to coordinate their evolution in the right 
way, even though they were never in causal contact? We must somehow modify the causal 
structure of the conventional FRW cosmology. 



5.3 Unwanted Relics 

Wc have spoken in the last lecture about grand unified theories (GUTs), and also about 
topological defects. If grand unification occurs with a simple gauge group G, any spontaneous 
breaking of G satisfies Tr2{G/H) — t:\{H) for any simple subgroup H. In particular, breaking 
down to the standard model will lead to magnetic monopoles, since 

T^2{G/H) = 7ri([SU(3) X SU(2) x \J{1)]/Ze) = Z . (186) 

(The gauge group of the standard model is, strictly speaking, [SU(3) x SU(2) x '[J{1)]/Zq. 
The Zq factor, with which you may not be familiar, only affects the global structure of the 
group and not the Lie algebra, and thus is usually ignored by particle physicists.) 

Using the Kibble mechanism, the expected relic abundance of monopoles works out to 

be 

This is far too big; the monopole abundance in GUTs is a serious problem for cosmology if 
GUTs have anything to do with reality. 

In addition to monopoles, there may be other model-dependent relics predicted by your 
favorite theory. If these arc incompatible with current limits, it is necessary to find some 
way to dilute their density in the early universe. 
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5.4 The General Idea of Inflation 



The horizon problem especially is an extremely serious problem for the standard cosmology 
because at its heart is simply causality. Any solution to this problem is therefore almost 
certain to require an important modification to how information can propagate in the early 
universe. Cosmological inflation is such a mechanism. 

Before getting into the details of inflation we will just sketch the general idea here. The 
fundamental idea is that the universe undergoes a period of accelerated expansion, defined 
as a period when a > 0, at early times. The effect of this acceleration is to quickly expand 
a small region of space to a huge size, diminishing spatial curvature in the process, making 
the universe extremely close to flat. In addition, the horizon size is greatly increased, so 
that distant points on the CMB actually are in causal contact and unwanted relics are 
tremendously diluted, solving the monopole problem. As an unexpected bonus, quantum 
fluctuations make it impossible for inflation to smooth out the universe with perfect precision, 
so there is a spectrum of remnant density perturbations; this spectrum turns out to be 
approximately scale-free, in good agreement with observations of our current universe. 

5.5 Slowly-Rolling Scalar Fields 

If inflation is to solve the problems of the standard cosmology, then it must be active at 
extremely early times. Thus, we would like to address the earliest times in the universe 
amenable to a classical description. We expect this to be at or around the Planck time tp 
and since Planckian quantities arise often in inflation we will retain values of the Planck mass 
in the equations of this section. There are many models of inflation, but because of time 
constraints we will concentrate almost exclusively on the chaotic inflation model of Linde. 
We have borrowed heavily in places here from the excellent text of Liddle and Lyth [7] . 

Consider modeling matter in the early universe by a real scalar fleld 0, with potential 
V{(f)). The energy- momentum tensor for (f) is 



For simplicity we will specialize to the homogeneous case, in which all quantities depend 
only on cosmological time t and set k — 0. A homogeneous real scalar fleld behaves as a 
perfect fluid with 



(188) 



\4>' - V{<P) . 



(189) 



(190) 



The equation of motion for the scalar fleld is given by 




(191) 
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which can be thought of as the usual equation of motion for a scalar field in Minkowski space, 
but with a friction term due to the expansion of the universe. The Friedmann equation with 
such a field as the sole energy source is 
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G 



+ v{<t>) 



(192) 



A very specific way in which accelerated expansion can occur is if the universe is dom- 
inated by an energy component that approximates a cosmological constant. In that case 
the associated expansion rate will be exponential, as we have already seen. Scalar fields can 
accomplish this in an interesting way. From (189) it is clear that if cjp' <^ 1^(0) then the 
potential energy of the scalar field is the dominant contribution to both the energy density 
and the pressure, and the resulting equation of state is p ~ — p, approximately that of a 
cosmological constant, the resulting expansion is certainly accelerating. In a loose sense, 
this negligible kinetic energy is equivalent to the fields slowly rolling down its potential; an 
approximation which we will now make more formal. 

Technically, the slow-roll approximation for inflation involves neglecting the (j) term 
in (191) and neglecting the kinetic energy of compared to the potential energy. The 
scalar field equation of motion and the Friedmann equation then become 

(193) 

//2 ^ , (194) 

where in this lecture a prime denotes a derivative with respect to 0. 

These conditions will hold if the two slow-roll conditions are satisfied. These are 

|e| < 1 

1^1 < 1 , (195) 



where the slow-roll parameters are given by 

Ml fV'V 



and 




(196) 



V = . (197) 

It is easy to see that the slow roll conditions yield infiation. Recall that infiation is defined 
by a/a > 0. We can write 

-^H + H\ (198) 
a 

so that infiation occurs if 

> -1 . (199) 



64 



But in slow-roll 



H 



(200) 



which will be small. Smallness of the other parameter rj helps to ensure that inflation will 
continue for a sufficient period. 

It is useful to have a general expression to describe how much inflation occurs, once it 
has begun. This is typically quantified by the number of e-folds, defined by 



Usually we are interested in how many efolds occur between a given field value (f) and the 
field value at the end of inflation 0cnd- deflned by e(0end) = 1- We also would like to express 
N in terms of the potential. Fortunately this is simple to do via 



The issue of initial conditions for inflation is one that is quite subtle and we will not get 
into a discussion of that here. Instead we will remain focused on chaotic inflation, in which 
we assume that the early universe emerges from the Planck epoch with the scalar fleld taking 
different values in different parts of the universe, with typically Planckian energies. There 
will then be some probability for inflation to begin in some places, and we shall focus on 
those. 

5.6 Attractor Solutions in Inflation 

For simplicity, let us consider a particularly simple potential 



where m has dimensions of mass. Wc shall also assume initial conditions that at the end of 




(201) 




(202) 




(203) 



the quantum epoch, which we label as t = 0, p ~ M^. 
The slow-roll conditions give that. Bit t — 0, 




(204) 



where we have deflned e = m/Mp <^ 1. We also have 



H 



(205) 



where 




(206) 
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so that the scalar field equation of motion becomes 

+ 3q;# + m^cf) = . (207) 

This is solved by 

= 00 - /9t , (208) 
where /3 = m^/Sa. Now, the slow-roll conditions remain satisfied provided that 

= (209) 



and therefore, (208) is vahd for a time 



At~^~^. (210) 



Why is this important? This is important because (208) is an attractor solution! Let's 
see how this arises. Consider perturbing (208) by writing 

(t>{t)^(t>{t) + x{t) , (211) 

where xit) is a "small" perturbation, substituting in to the equation of motion and linearizing 
in X- One obtains 

X + 3a0x = . (212) 

This equation exhibits two solutions. The first is just a constant, which just serves to 
move along the trajectory. The second solution is a decaying mode with time constant 
tc — (3a0)~^. Since tc -C At, all solutions rapidly decay to (208) - it is an attractor. 

5.7 Solving the Problems of the Standard Cosmology 

Let's stick with the simple model of the last section. It is rather easy to see that the Einstein 
equations are solved by 



a{t) = a(0) exp 



p 



(213) 



as we might expect in inflation. The period of time during which (208) is valid ends at 
t At, at which 

a(At) ~ a(0) exp (^^) , (214) 

Now, taking a typical value for m, for which e < 10^'^, we obtain 

a{At) ~ a(0) X lO^-^xio* . (215) 

This has a remarkable consequence. A proper distance Lp at t — will inflate to a size 
10^° cm after a time At ~ 5 x 10~^^ seconds. It is important at this point to know that 
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the size of the observable universe today is Hq^ ~ andhasnothadenoug hlO^^cml Therefore, 
only a small fraction of the original Planck length comprises today's entire visible universe. 
Thus, homogeneity over a patch less than or of order the Planck length at the onset of 
inflation is all that is required to solve the horizon problem. Of course, if we wait sufficiently 
long we will start to see those inhomogeneities (originally sub-PIanckian) that were inflated 
away. However, if inflation lasts long enough (typically about sixty e-folds or so) then this 
would not be apparent today. Similarly, any unwanted rehcs are diluted by the tremendous 
expansion; so long as the GUT phase transition happens before inflation, monopoles will 
have an extremely low density. 

Inflation also solves the flatness problem. There are a couple of ways to see this. The 
first is to note that, assuming inflation begins, the curvature term in the Friedmann equation 
very quickly becomes irrelevant, since it scales as a{t)~^. Of course, after inflation, when 
the universe is full of radiation (and later dust) the curvature term redshifts more slowly 
than both these components and will eventually become more important than both of them. 
However, if inflation lasts sufficiently long then even today the total energy density will be 
so close to unity that wc will notice no curvature. 

A second way to see this is that, following a similar analysis to that leading to (173) 
leads to the conclusion that Q = 1 is an attractor rather than a repeller, when the universe 
is dominated by energy with w < —1/3. Therefore, Jl is forced to be very close to one by 
inflation; if the inflationary period lasts sufficiently long, the density parameter will not have 
had time since to stray appreciably away from unity. 



5.8 Vacuum Fluctuations and Perturbations 

Recall that the structures - clusters and superclusters of galaxies - we see on the largest 
scales in the universe today, and hence the observed fluctuations in the CMB, form from the 
gravitational instability of initial perturbations in the matter density. The origin of these 
initial fluctuations is an important question of modern cosmology. 

Inflation provides us with a fascinating solution to this problem - in a nutshell, quantum 
fluctuations in the inflaton flcld during the inflationary epoch are stretched by inflation and 
ultimately become classical fluctuations. Let's sketch how this works. 

Since inflation dilutes away all matter flelds, soon after its onset the universe is in a 
pure vacuum state. If we simplify to the case of exponential inflation, this vacuum state is 
described by the Gibbons-Hawking temperature 

H VV , , 

where we have used the Friedmann equation. Because of this temperature, the inflaton 
experiences fluctuations that are the same for each wavelength = Tgh- Now, these 
fluctuations can be related to those in the density by 

^p=^.0. (217) 
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Inflation therefore produces density perturbations on every scale. The amphtude of the 
perturbations is nearly equal at each wavenumber, but there will be slight deviations due to 
the gradual change in V as the inflaton rolls. We can characterize the fluctuations in terms 
of their spectrum As{k), related to the potential via 



Alik) ^ 



(218) 

k=aH 



where k = aH indicates that the quantity V^/{V')'^ is to be evaluated at the moment when 
the physical scale of the perturbation A = a/A; is equal to the Hubble radius H~^. Note 
that the actual normalization of (218) is convention-dependent, and should drop out of any 
physical answer. 

The spectrum is given the subscript "S" because it describes scalar fluctuations in the 
metric. These are tied to the energy-momentum distribution, and the density fluctuations 
produced by inflation are adiabatic — fluctuations in the density of all species are corre- 
lated. The fluctuations are also Gaussian, in the sense that the phases of the Fourier modes 
describing fluctuations at different scales are uncorrelated. These aspects of inflationary per- 
turbations — a nearly scale-free spectrum of adiabatic density fluctuations with a Gaussian 
distribution — are all consistent with current observations of the CMB and large-scale struc- 
ture, and have been conflrmed to new precision by WMAP and other CMB measurements. 

It is not only the nearly- massless inflaton that is excited during inflation, but any ncarly- 
massless particle. The other important example is the graviton, which corresponds to tensor 
perturbations in the metric (propagating excitations of the gravitational field) . Tensor fiuc- 
tuations have a spectrum 



Aifk) ~ 



(219) 

k=aH 

The existence of tensor perturbations is a crucial prediction of inflation which may in prin- 
ciple be veriflable through observations of the polarization of the CMB. Although CMB 
polarization has already been detected [166], this is only the ii^-mode polarization induced 
by density perturbations; the B-mode polarization induced by gravitational waves is ex- 
pected to be at a much lower level, and represents a signiflcant observational challenge for 
the years to come. 

For purposes of understanding observations, it is useful to parameterize the perturbation 
spectra in terms of observable quantities. We therefore write 

Al{k) oc (220) 

and 

Al{k) (X A;"^ , (221) 

where ns and nx are the "spectral indices" . They are related to the slow-roll parameters of 
the potential by 

ns = 1 - 6e + 277 (222) 
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and 

riT = -2e . (223) 

Since the spectral indices are in principle observable, we can hope through relations such as 
these to glean some information about the inflaton potential itself. 

Our current knowledge of the amplitude of the perturbations already gives us important 
information about the energy scale of inflation. Note that the tensor perturbations depend 
on V alone (not its derivatives), so observations of tensor modes yields direct knowledge 
of the energy scale. If large-scale CMB anisotropics have an appreciable tensor component 
(possible, although unlikely), we can instantly derive Vinflation ~ (10^^ GeV)^. (Here, the 
value of V being constrained is that which was responsible for creating the observed fluc- 
tuations; namely, 60 e- folds before the end of inflation.) This is remarkably reminiscent of 
the grand unification scale, which is very encouraging. Even in the more likely case that the 
perturbations observed in the CMB are scalar in nature, we can still write 

^n£tion - e'^Ho^' GeV , (224) 

where e is the slow-roll parameter defined in (196). Although we expect e to be small, the 
1/4 in the exponent means that the dependence on e is quite weak; unless this parameter is 
extraordinarily tiny, it is very likely that l^„flttion ~ 10^^-10^^ GeV. 

We should note that this discussion has been phrased in terms of the simplest models of 
inflation, featuring a single canonical, slowly-rolling scalar field. A number of more complex 
models have been suggested, allowing for departures from the relations between the slow-roll 
parameters and observable quantities; some of these include hybrid infiation [167, 168, 169], 
infiation with novel kinetic terms [170]. the curvaton model [171, 172, 173], low-scale models 
[174, 175], brane infiation [176, 177, 178, 179, 180, 181, 182] and models where perturbations 
arise from modulated coupling constants [183, 184, 185, 186, 187]. This list is necessarily 
incomplete, and continued exploration of the varieties of inflationary cosmology will be a 
major theme of theoretical cosmology into the foreseeable future. 

5.9 Reheating and Preheating 

Clearly, one of the great strengths of inflation is its ability to redshift away all unwanted 
relics, such as topological defects. However, inflation is not discerning, and in doing so any 
trace of radiation or dust-like matter is similarly redshifted away to nothing. Thus, at the 
end of infiation the universe contains nothing but the infiationary scalar field condensate. 
How then does that matter of which we are made arise? How does the hot big bang phase 
of the universe commence? How is the universe reheated^. 

For a decade or so reheating was thought to be a relatively well-understood phenomenon. 
However, over the last 5-10 years, it has been realized that the story may be quite compli- 
cated. We will not have time to go into details here. Rather, we will sketch the standard 
picture and mention briefiy the new developments. 

Inflation ends when the slow-roll conditions are violated and, in most models, the fleld 
begins to fall towards the minimum of its potential. Initially, all energy density is in the 
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inflaton, but this is now damped by two possible terms. First, the expansion of the universe 
naturally damps the energy density. More importantly, the inflaton may decay into other 
particles, such as radiation or massive particles, both fermionic and bosonic [188, 189]. To 
take account of this one introduces a phenomenological decay term into the scalar field 
equation. If we focus on fermions only, then a rough expression for how the energy density 
evolves is 



The inflaton undergoes damped oscillations and decays into radiation which equilibrates 
rapidly at a temperature known as the reheat temperature Trh- 

In the case of bosons however, this ignores the fact that the inflaton oscillations may 
give rise to parametric resonance. This is signified by an extremely rapid decay, yielding a 
distribution of products that is far from equilibrium, and only much later settles down to an 
equilibrium distribution at energy Trh- Such a rapid decay due to parametric resonance is 
known as preheating [190, 191]. 

One interesting outcome of preheating is that one can produce particles with energy far 
above the ultimate reheat temperature. Thus, one must beware of producing objects, such as 
topological relics, that inflation had rid us of. However, preheating can provide some useful 
effects, such as interesting new ways to generate dark matter candidates and topological 
relics and to generate the baryon asymmetry [192, 193, 194, 195, 196, 197]. 

5.10 The Beginnings of Inflation 

If inflation begins, we have seen that it achieves remarkable results. It can be shown quite 
generally that the start of inflation in a previously non-inflating region requires pre-existing 
homogeneity over a horizon volume [198]. It seems unlikely that classical, causal physics 
could cause this and we may interpret this as an indication that quantum physics may play 
a role in the process of inflation. 

Although it can be shown [199, 200, 201, 202, 203] that inflation cannot be past eternal, 
one consequence of quantum mechanics is that inflation may be future eternal [204, 205, 
206, 207, 208] . To see this, consider inflation driven by a scalar fleld rolling down a simple 
potential as in flgure (5.10). 

Imagine that, in a given space-like patch, inflation is occurring with the scalar fleld 
at a potential energy Vi. Although the classical motion of the scalar fleld is towards the 
minimum of the potential, there is a nonzero probability that quantum fluctuations will 
cause the inflaton to jump up to a higher value of the potential V2 > Vi over a sufficiently 
large part of the inflating patch. In the new region, the Fricdmann equation tells us that 
the local Hubble parameter is now larger than in the inflating background 



and requiring the region to be sufficiently large means that its linear size is at least H2 , 
which may be considerably smaller than the size of the background inflationary horizon. 



+ {3H + T^)p^ = . 



(225) 



H2 




(226) 
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Thus, inflation occurs in our new region with a larger Hubble parameter and, because of the 
exponential nature of inflation, the new patch very quickly contains a much larger volume 
of space than the old one. 

There are, of course, parts of space in which the classical motion of the infiaton to its 
minimum completes, and those regions drop out of the inflationary expansion and are free 
to become regions such of our own. However, it is interesting to note that, as long as the 
probability for inflation to begin at all is not precisely zero, any region like ours observed 
at a later time comes from a previously inflating region with probability essentially unity. 
Nevertheless, in this picture, most regions of the universe inflate forever and therefore this 
process is known as eternal inflation. 

If correct, eternal inflation provides us with an entirely new global view of the universe and 
erases many of our worries about the initial condition problems of inflation. However, some 
caution is in order when considering this model. Eternal inflation relies on the quantum 
backreaction of geometry to matter fluctuations over an extended spatial region. This is 
something we do not know how to treat fully at the present time. Perhaps string theory can 
help us with this. 
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